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1  Introduction 


This  report  describes  the  background  theory,  experimental  techniques,  and  results 
of  research  relating  to  iidar  polarization,  scattering,  and  absorption  effects  which 
PhotoMetrics  performed,  under  Contract  F19628-86-C-0102,  using  the  GL/OP  low 
altitude  lidar  system.  A  companion  report,  hereafter  in  this  report  referred  to  as 
Part  II,  describes  the  upgrades  and  modifications  PhotoMetrics  made  and  provides 
an  instruction  manual  for  the  operation  of  the  lidar  system. 

PhotoMetrics  has  used  the  low  altitude  lidar  system  to  gather  data  in  support 
of  GL’s  broad  program  aimed  at  modeling  the  optical  properties  of  the  lower  at¬ 
mosphere.  The  general  approach  we  have  taken  in  our  research  efforts  has  been  in 
improving  the  understanding  of  lidar  as  a  diagnostic  tool  for  studying  the  atmo¬ 
sphere.  We  have  concentrated  this  approach  in  three  broad  areas: 

1.  understanding  the  effects  of  multiple  scattering  in  lidar  experiments, 

2.  understanding  the  effects  of  scattering  from  irregularly  shaped,  randomly  ori¬ 
ented  particles,  and 

3.  determining  the  validity  and  improving  upon  techniques  of  calculating  the 
extinction  coefficient  from  lidar  data. 

Sections  3  to  6  each  report  on  major  tasks  undertaken  during  this  program.  Each 
has  a  bearing  on  one  of  the  above  three  mentioned  areas. 

Section  2  gives  a  brief  description  of  the  lidar  system  and  its  measurement 
capabilities.  A  detailed  description  is  given  in  Part  II. 

Section  3  describes  our  measurements  of  the  polarization  state  of  the  backscat- 
tered  532  nm  radiation.  Most  of  these  measurements  have  been  done  during  falling 
snow  conditions. 

Section  4  gives  a  brief  description  of  Mie  scattering  and  discusses  the  results  of 
a  Mie  scattering  program  that  we  have  written  for  an  MS-DOS  based  computer. 
This  program  was  used  in  conjunction  with  understanding  the  validity  of  the  lidar 
calculated  extinction  coefficient. 
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Section  5  describes  our  work  on  extinction.  It  describes  the  standard  extinction 
algorithm,  shows  results  when  applied  to  lidar  data,  and  describes  some  enhance¬ 
ments  we  have  made. 

Section  6  discusses  our  preliminary  efforts  of  modeling  multiple  scattering  effects 
in  a  lidar  experiment. 

Appendix  A  discusses  and  shows  results  of  the  field  trip  to  the  desert  region  of 
Dripping  Spring,  New  Mexico  in  June,  1988.  This  trip  was  made  in  support  of  a 
measurement  program  at  Sandia  National  Laboratory.  This  section  also  discusses 
the  minimum  calculable  extinction  coefficient  as  a  function  of  the  characteristics  of 
a  lidar  system.  These  calculations  were  motivated  by  the  Dripping  Spring  field  trip. 

2  Lidar  Parameters 

The  GL/OP  mobile,  low  altitude  lidar  is  a  biaxial  system  with  the  receiver  mounted 
adjacent  to  the  transmitter.  The  combination  is  housed  in  a  temperature  controlled 
enclosure  which  is  mounted  on  a  steerable  trunion  on  top  of  a  trailer.  The  data 
acquisition  and  control  systems  are  located  inside  the  trailer. 

The  laser  transmitter  is  a  Nd:YAG  operating  at  20  Hz.  The  primary  output  at 
1.06  nm  is  120  mJ/pulse.  With  a  temperature-tuned  CD*A  crystal  25  mJ/pulse 
of  the  primary  is  converted  to  0.53  nm.  The  laser  pulse  is  approximately  15  ns 
in  length.  The  beam  divergence  is  approximately  2.5  mrad.  The  output  is  nearly 
100%  polarized. 

The  receiver  is  an  f/ 1.3  refractive  telescope  with  a  15  cm  diameter  aperture  and 
a  10  mrad  field  of  view.  The  transmit, ter /receiver  crossover  point  is  at  150  m.  A 
gated  photomultiplier  tube  detects  0.53  /xm  radiation  through  a  10A  bandwidth 
filter.  For  most  of  the  work  described  in  this  report  the  PMT  was  operated  in 
current  mode,  although  it  will  simultaneously  operate  in  photon  counting  mode 
(described  in  Part  II).  A  silicon  avalanche  photodiode  detects  1.06  /xm  radiation 
through  a  10  nm  bandwidth  filter.  The  polarization  state  of  the  backscattered  0.53 
H m  radiation  can  be  determined  with  a  polarizer  wheel  located  inside  the  receiver. 

The  data  acquisition  system  consists  of  a  CAMAC  crate  interfaced  to  an  80386 
computer.  For  each  channel,  signal  averaging  is  performed  by  a  special  CAMAC  av- 
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eraging  memory  directly  interfaced  to  a  transient  recorder.  The  transient  recorders 
are  12  bit  with  an  8192  word  memory  and  can  acquire  data  at  rates  up  to  20  Mhz. 
We  typically  run  the  transient  recorders  at  10  MHz.  The  signal  averagers  are  24 
bit  and  will  accumulate  data  for  up  to  65536  laser  shots. 

3  Polarization  Studies 

3.1  Introduction 

PhotoMetrics  designed  and  implemented  a  system  to  measure  the  polarization  state 
of  the  backscattered  532  nm  radiation.  The  details  of  the  system  and  its  method  of 
operation  are  given  in  Part  II. 

Here  we  describe  the  polarization  studies  that  we  performed.  Our  most  exten¬ 
sive  work  was  done  for  falling  snow.  Section  3.2  provides  the  theoretical  background 
for  the  polarization  measurements.  Section  3.3  describes  the  general  qualitative  fea¬ 
tures  of  the  polarization  state  of  backscattered  lidar  radiation.  Section  3.4  describes 
our  studies  of  falling  snow,  both  experimental  and  theoretical. 

Knowledge  of  the  polarization  state  of  the  radiation,  together  with  multiple 
scattering  models,  is  necessary  in  order  to  determine  the  true  extinction  coeffi¬ 
cient.  A  medium  in  which  there  is  a  large  amount  of  multiple  scattering  has  a 
smaller  extinction  coefficient  than  would  otherwise  be  calculated  by  assuming  sin¬ 
gle  scattering  (there  is  more  backscattered  radiation  measured  than  is  accounted  for 
theoretically).  Furthermore,  since  there  are  more  multiple  scatters  as  the  radiation 
penetrates  deeper  into  the  medium,  the  extinction  calculation  using  single  scatter 
assumptions  becomes  less  accurate. 

Knowledge  of  the  degree  of  polarization,  together  with  previous  experiments  or 
models,  can  be  used  to  identify  the  scattering  particles.  The  degree  of  polarization 
of  backscattered  radiation  from  a  medium  composed  of  irregularly  shaped  particles 
is  characteristic  of  the  particular  particles.  For  example,  snowflakes  will  induce  a 
greater  amount  of  depolarization  than  the  hexagonally  shaped  ice  crystals  of  cirrus 
clouds  because  of  the  wider  distribution  of  randomly  oriented  surfaces  within  a 
snowflake.  The  radiation  in  snow  generally  reflects  many  more  times  off  several 
surfaces  of  the  same  snowflake  before  scattering.  Thus,  single  scattering  from  a 
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snowflake  can  be  viewed  as  multiple  scattering  within  that  snowflake. 

3.2  Theoretical  Background 

The  polarization  state  of  electromagnetic  radiation  (coherent  or  incoherent,  monochro¬ 
matic  or  nonmonchromatic)  is  completely  specified  by  its  four  Stokes  parameters. 
These  four  parameters  (the  set  is  not  unique)  describe 

l)  the  fraction  of  radiation  which  is  in  a  single  polarization  state  (i.e., 
the  degree  of  polarization), 

2,3)  the  amplitudes  of  two  perpendicular  components  of  the  electric  (or 
magnetic)  field  of  the  polarized  portion  of  the  radiation,  and 
4)  the  phase  between  these  two  components. 

For  a  quasi-monochromatic  wave  described  by 

i  =  [F;I(t)x  +  F;v(0j/]e,(“-w‘)  (i) 

where 

Et  =  «,(t)e*'W 

Ev  =  ev(t)e»>(‘)  (2) 

and  Ex  v(t)  are  slowly  varying  compared  to  e~lut,  four  commonly  used  Stokes  pa¬ 
rameters  are  [l] 

I  =  (EtE:)  +  {EyE;)  =  (el  +  el) 

Q  =  (EtE:)-(EvE'v)  =  (el-el) 

U  =  {EZE'V)  +  (E*tEy)  =  (2exevcos(0I  -  <pv )) 

V  =  *  [(EtE'v)  -  <£;£„)]  =  (2cIe„sin(<^I  -  <£„)>.  (3) 

The  symbol  (•  •  •)  denotes  time  average  and  denotes  complex  conjugation.  In 
addition,  it  was  assumed  that  ex  y  and  <f>XtV  are  real  quantities.  The  intensities  of 
light  polarized  in  the  x  and  y  directions  are  6*  and  £*,  respectively. 

If  the  radiation  is  completely  polarized  then 

Q2  +  U2  +  V2  =  1 2  (4) 
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and  there  are  therefore  only  three  independent  parameters.  More  generally,  a  por¬ 
tion  of  the  radiation  has  random  polarization  which  gives  the  relation 

Q2  +U2  +  V2  <  I2.  (5) 


A  useful  parameter  is  the  degree  of  polarization  defined  as 

y/Q2  +  U*  +  V2  Ip 

I  /p  +  /«’ 


(6) 


where  Ip  —  y/Q 2  +  U2  +  V2  is  the  intensity  of  the  polarized  portion  and  /„  =  I  —  Ip 
is  the  intensity  of  the  unpolarized  portion  of  the  radiation. 

It  is  illuminating  to  visualize  the  normalized  quantities  Q/I,  U/I,  and  V /I  as 
axes  of  a  three  dimensional  cartesian  coordinate  system,  as  shown  in  Fig.  1.  The 
polarization  state  of  the  radiation  is  then  represented  by  a  vector  (f ,  7)  which 

lies  inside  or  on  the  unit  sphere  (referred  to  as  the  Poincare  sphere  [2]).  The  length 
of  the  vector  is  the  degree  of  polarization.  The  component  of  the  vector  in  the 
^  —  j  plane  is  the  degree  of  linear  polarization  (defined  as  \Q 2  +  U2f2/I).  The 
component  of  the  vector  along  the  y  axis  is  the  degree  of  circular  polarization. 
Other  quantities,  such  as  ellipticity  and  azimuth  of  the  polarized  portion  of  the 
radiation  are  also  easily  visualized. 

An  important  property  of  the  Stokes  parameters  for  lidar  work  is  that  the  Stokes 
parameters  of  the  total  radiation  equal  the  sum  of  the  Stokes  parameters  for  the 
individual  parts.  With  lidar  we  observe,  during  each  digitization  period  of  the  data 
acquisition  system,  the  total  backscattered  radiation  from  many  scatters  within 
a  volume  of  atmosphere.  The  light  from  each  scatterer  is  described  by  its  own 
set  of  Stokes  parameters.  At  the  receiver  all  the  backscattered  radiation  combines 
and  has  a  new  set  of  Stokes  parameters.  If  radiation  from  the  different  scatterers 
is  uncorrelated,  as  is  usually  the  case  for  atmospheric  scattering,  then  the  Stokes 
parameters  are  additive: 


l 1 ) 

(  i 

Q 

u 

=  £ 

Qi 

ut 

U  J 

{  v 

The  set  of  four  measurements  required  to  determine  the  four  Stokes  parameters 
of  the  backscattered  radiation  is  not  unique.  In  the  low  altitude  lidar  system  we 
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circular 
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Figure  1:  The  Poincar^  sphere:  a  geometrical  interpretation  of  the  Stokes  parame¬ 
ters. 
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make  the  following  four  measurements  (referenced  to  the  laser  output,  which  is 
nearly  100%  linearly  polarized  and  represented  by  the  Stokes  vector  [7,  7,0, Oj): 

7)1 :  intensity  polarized  parallel  to  the  laser  output, 

I±:  intensity  polarized  perpendicular  to  the  laser  output, 

/is:  intensity  polarized  45°  to  laser  output,  and 
7q:  one  of  the  circularly  polarized  components. 

The  Stokes  parameters  of  the  backscattered  light  are  then 

7  =  7||  +  7X 

Q  =  h-h 

U  =  2746  —  7),  —  7_l 

V  =  270  —  7||  —  7X  (8) 

A  scattering  event  is  described  as  an  incident  4x1  Stokes  vector  being  trans¬ 
formed  into  a  scattered  4x1  Stokes  vector.  If  the  scattering  is  linear  (i.e.,  if  the 
scattering  does  not  depend  on  the  electric  or  magnetic  field  strengths,  which  is  gen¬ 
erally  the  case  for  situations  of  interest  to  us)  then  a  4  x  4  matrix  (the  scattering 
matrix),  which  is  a  function  of  scattering  angle,  represents  the  transformation.  Only 
seven  independent  parameters  are  needed  to  specify  the  sixteen  component  matrix. 

3.3  General  Behavior  of  Backscattered  Radiation 

Backscattered  radiation  from  spherical  particles  of  all  sizes  (e.g.,  many  types  of 
aerosols  and  water  droplets  in  cumulus  clouds)  and  Rayleigh  backscattered  radia¬ 
tion  (i.e.,  scattering  from  particles  small  compared  to  the  wavelength)  retain  the 
polarization  of  the  incident  light.1  In  such  cases  the  cylindrical  symmetry  about  the 
axis  defined  by  the  laser  beam  prevents  coupling  to  a  different  polarization  state. 

The  backscattered  radiation  received  by  a  lidar  changes  its  polarization  state 
relative  to  that  of  the  incident  radiation  if 

•  the  radiation  multiple  scatters  or 

1  For  Rayleigh  scattering  from  particles  with  nonsero  dipole  moment  the  polarisation  state  of 
backscattered  radiation  can  change  relative  to  the  incident  radiation.  However,  this  is  usually  a 
small  effect  and  is  neglected  in  lidar  work. 
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•  the  radiation  single  scatters  from  many  irregularly  shaped,  randomly  oriented 
particles. 

The  first  situation  occurs  in  most  clouds,  where  the  optical  thickness  (extinction 
coefficient  times  path  length)  is  greater  than  one.  The  second  situation  occurs  in 
cirrus  clouds  (hexagonal  ice  crystals),  falling  snow,  and  with  some  types  of  aerosols. 

In  both  situations,  if  the  laser  output  is  completely  polarized  then  the  backscat- 
tered  radiation  generally  is  partially  polarized,  or  has  a  degree  of  polarization  less 
than  one.  The  polarized  portion  of  the  radiation  will,  on  average  (i.e.,  for  many 
laser  pulses),  be  in  the  same  polarization  state  as  the  incident  light  (with  some 
spread  about  this  average).  From  the  point  of  view  of  the  Poincare  sphere,  the 
scattered  vector  is  parallel  to  the  incident  vector,  but  smaller  in  magnitude. 

For  these  situations  two  measurements,  instead  of  four  measurements,  are  suffi¬ 
cient  to  to  determine  the  polarization  state  of  the  backscattered  light.2  For  linearly 
polarized  incident  light,  represented  by  vector  (1,0,0)  on  the  Poincare  sphere,  the 
scattered  vector  would  be  {Q/1, 0,0),  where  Q/I  <  1.  The  quantity  Q/I  in  this  sit¬ 
uation  is  both  the  degree  of  polarization  and  the  degree  of  linear  polarization.  Thus, 
two  sufficient  measurements  may  be  the  polarization  components  of  the  backscat¬ 
tered  light  parallel  and  perpendicular  to  the  polarization  of  the  incident  light.  The 
values  I  and  Q  can  be  determined  from  these  two  measurements. 

The  spread  in  the  polarization  state  about  the  average  value,  mentioned  above, 
may  be  a  valuable  piece  of  information  needed  to  characterize  the  scattering  medium. 
For  scattering  from  irregularly  shaped,  randomly  oriented  particles,  the  spread  rep¬ 
resents  the  fact  that  each  single  scatter  may  change  the  polarization  state  slightly, 
but  the  random  orientation  does  not  allow  a  preferred  shift  of  the  average.  The 
size  of  the  spread  is  reduced  in  proportion  to  the  square  root  of  (the  number  of 
scatterers  in  the  receiver  viewing  volume  multiplied  by  the  number  of  laser  shots 
in  the  lidar  run).  The  complete  set  of  four  polarization  measurements  is  needed  to 
determine  the  spread.  We  have  not  made  use  of  this  information  in  our  studies. 

2 Four  measurements  are  only  needed  if  the  polarized  portion  of  the  scattered  radiation  is  in 
a  different  state  from  that  of  the  incident  radiation,  such  as  for  single  scattering  from  a  single 
irregularly  shaped  particle,  or  from  many  irregularly  shaped,  oriented  particles. 
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3.4  Falling  Snow 

3.4.1  Motivation 

The  identification  of  cirrus  clouds,  particularly  high,  thin,  subvisible  cirrus  clouds, 
was  the  impetus  for  implementing  the  polarization  detection  system.  As  described 
in  Section  3.3,  the  backscattered  radiation  from  the  hexagonal  shaped  ice  crystals 
of  cirrus  is  expected  to  have  an  unpolarized  component.  Further  depolarization 
occurs  due  to  multiple  scattering. 

Unfortunately,  clouds  are  relatively  quickly  varying  media.  We  have  observed 
that  many  types  of  clouds  will  either  change  their  physical  shapes  or  drift  within  a 
few  seconds.  A  full  polarization  measurement,  consisting  of  four  5-10  second  lidar 
runs  (100-200  laser  shots),  usually  cannot  be  made  within  this  time.  Even  a  two 
lidar  run  polarization  measurement  takes  too  long,  unless  we  are  willing  to  sacrifice 
signal-to-noise  by  averaging  over  fewer  laser  shots.  The  best  we  can  hope  for,  in  the 
current  setup,  is  to  get  an  average  polarization  measurement  through  the  cloud.8 

Falling  snow  is  an  excellent  alternative  to  clouds  in  studying  polarization  ef¬ 
fects.  Falling  snow  is  a  relatively  slowly  varying  medium.  Snow  density  may  vary 
as  quickly  as  a  few  seconds,  but  as  slowly  as  several  minutes.  Falling  snow  typically 
has  a  smaller  extinction  coefficient  than  clouds,  providing  a  longer  range  over  which 
single  scattering  can  be  observed  separately  from  multiple  scattering  effects.  The  ir¬ 
regularly  shaped,  randomly  oriented  snowflakes  display  the  depolarization  effects  of 
single  scattered  light  quite  nicely.  In  addition,  since  snow  falls  to  the  ground,  other 
ground  based  diagnostics  may  be  used  to  collect  companion  data  for  correlation  and 
additional  information. 

Sections  3.4.2  to  3.4.5  describe  our  experimental  and  theoretical  work  on  snow 
scattering.  Section  3.4.2  presents  observations  and  conclusions  from  lidar  data 
acquired  during  falling  snow  conditions.  Section  3.4.3  discusses  our  work  on  deter¬ 
mining  the  maximum  range  of  usefulness  of  lidar  data  for  polarization  calculations. 
Section  3.4.4  presents  a  scenario  for  snowflake  scattering  that  we  have  formulated 
based  our  observations  and  conclusions.  Section  3.4.5  describes  the  modeling  that 

3  We  have  designed  a  polarization  system,  to  be  implemented  in  a  new  lidar  receiver,  which  makes 
near  simulaneous  measurements  of  the  four  polarisations.  This  system  is  described  in  Part  II. 
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we  have  done  to  describe  single  snowflake  scattering. 

3.4.2  Experiment 

Lidar  experiments  performed  in  falling  snow  indicate  that 

•  the  degree  of  polarization  of  backscattered  radiation  is  in  the  range  of  0.3-0.5 
and  is  fairly  insensitive  to  snow  density  and 

•  the  polarized  portion  of  the  backscattered  radiation  is,  on  average,  in  the  same 
polarization  state  as  the  incident  radiation. 

Below  we  present  the  results  of  data  analysis  which  support  these  conclusions.  We 
also  present  a  scenario  for  snowflake  scattering  based  on  these  data  and  conclusions. 

Figure  2a  shows  four  successive  lidar  returns  from  falling  snow.  Each  return 
corresponds  to  one  of  the  four  polarization  measurements  of  Eq.  8. 

These  snow  returns  are  four  of  several  hundred  returns  taken  over  a  four  hour 
period  on  Feb.  4,  1988.  During  this  period  the  snowfall  rate,  calculated  from  data 
obtained  with  a  Particle  Measurement  System  (PMS)  rain  distrometer,  varied  from 
zero  to  80  cm-Jsec-1, 

Figure  2b  shows  the  Stokes  parameters  and  degree  of  polarization  calculated 
from  these  data.  For  all  data  acquired  on  Feb  4,  1988  these  calculations  are  per¬ 
formed  in  the  manner  prescribed  by  Eqs.  8  and  6. 

For  comparison  Fig.  3  shows  the  Stokes  parameters  calculated  from  data  ac¬ 
quired  during  clear  air  conditions.  The  degree  of  polarization  for  clear  air  is  near 
unity  because  most  of  the  scatterers  are  spherically  shaped  aerosols. 

Figures  4  through  6  show  plots  of  the  Stokes  parameters  and  the  degree  of 
polarization  versus  snow  density  and  versus  extinction  coefficient  at  different  ranges 
for  all  the  lidar  returns  of  Feb.  4,  1988.  Snow  density  is  calculated  from  snow  rate 
obtained  from  the  rain  distrometer  by  assuming  a  snow  velocity  of  50  cm/sec. 
The  extinction  coefficient  is  determined  using  the  Klett  inversion  routine  on  the 
lidar  data.  We  discuss  these  calcualations,  and  their  validity,  in  more  detail  in 
Section  5.3.2  on  page  49. 
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Figure  2:  a)  Four  successive  lidar  returns  from  falling  snow,  b)  Stokes  parameters 
and  degree  of  polarization  calculated  from  these  data. 
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Stokes  parameters  for  "clear"  air 


Figure  3:  Stokes  parameters  and  the  degree  of  polarization  calculated  from  data 
during  clear  air  conditions. 
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Figure  7  shows  a  plot  of  the  range-averaged  value  of  Q/I  versus  the  range- 
average  value  of  extinction.  The  averages  are  computed  over  the  interval  200-500 
m.  The  value  Q/I  is  equal  to  the  degree  of  linear  polarization  when  the  polarized 
portion  of  the  radiation  is  in  the  same  polarization  state  as  that  of  the  incident 
laser  radiation. 

Figure  8  shows  plots  of  the  parameter  Q/I  versus  lidar  run  for  several  different 
ranges.  All  of  the  lidar  runs  for  Feb.  4,  1988  are  included. 

Figures  4  to  7  show  that  the  degree  of  polarization  is  insensitive  to  snow  density 
or,  equivalently,  extinction  coefficient.  Furthermore,  the  parameters  U //  and  V/I 
do  not  contribute  much  to  the  degree  of  polarization,  since  the  value  Q/I  is  nearly 
equal  to  the  degree  of  polarization.  Therefore,  the  degree  of  linear  polarization  is 
approximately  equal  to  the  degree  of  polarization,  which  is  approximately  equal  to 
Q/I.  (Put  another  way,  the  degree  of  circular  polarization  is  small). 

Since  the  incident  laser  radiation  is  entirely  Q/I  (i.e.,  U/I  =  V/I  =  0)  then 
the  backscattered  radiation  is  merely  depolarized,  or  given  a  random  component 
of  polarization.  The  part  that  remains  polarized  is  in  the  same  state  as  the  in¬ 
cident  radiation.  This  is  true  on  average,  since  there  is  some  spread  in  degree  of 
polarization  and  Q/I. 

Since  the  parameter  Q/I  is  approximately  equal  to  the  degree  of  polarization  we 
can  use  it  as  an  indicator  of  the  effects  of  snow  on  polarization  of  the  backscattered 
radiation.  Therefore,  when  discussing  snow,  we  henceforth  use  the  term  degree  of 
polarization  to  refer  to  the  value  Q/I. 

The  insensitivity  of  degree  of  polarization  is  also  revealed  by  data  acquired 
during  the  falling  snow  of  Jan  29,  1990.  Figure  9a  is  a  plot  of  range-averaged  Q/I 
versus  range-averaged  extinction  for  that  day.  The  snow  fall  on  this  day  varied  from 
light  to  heavy  (compare  extinction  coefficients  of  Jan.  29,  1990  to  those  of  Feb.  4, 
1988.) 

Figure  9b  is  the  same  plot  but  for  data  acquired  during  the  very  light  snow  of  Jan 
15,  1990.  The  averages  in  Figs.  9a  and  9b  are  performed  over  the  interval  200-400 
m.  The  extinction  coefficients  are  determined  with  the  Klett  inversion  routine. 

The  Jan.  29,  1990  data  shows  that  the  degree  of  polarization  is  insensitive 
to  extinction  coefficient  for  extinctions  greater  than  approximately  0.5-1  km-1. 
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a  (1/km)  (Klett) 


Figure  4:  Stokes  parameters  and  degree  of  polarization  versus  a)  snow  density  and 
b)  extinction  coefficient  at  150  m  range. 
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Figure  5:  Stokes  parameters  and  degree  of  polarization  versus  a)  snow  density  and 
b)  extinction  coefficient  at  285  m  range. 
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Figure  6:  Stokes  parameters  and  degree  of  polarization  versus  a)  snow  density  and 
b)  extinction  coefficient  at  435  m  range. 
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Figure  7:  Average  Q/I  v i 
4,  1988  data. 
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Figure  9:  Average  Q/I  versus  average  extinction  for  a)  Jan.  29,  1990  (light  to 
heavy  snow  and  b)  Jan.  15,  1990  (light  snow). 
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The  Jan.  15,  1990  data  does  not  have  extinctions  greater  than  these  values  and, 
therefore,  it  does  not  show  the  insensitivity  of  degree  of  polarization. 

This  sensitivity  of  degree  of  polarization  at  low  snow  density  is  explained  by  the 
fact  that  there  are  other  constituents  in  the  atmosphere  during  falling  snow  that 
do  not  depolarize  the  backscattered  radiation.  As  the  number  of  snow  scatterers 
decrease,  the  fraction  of  these  coherent  scatterers  increase.  We  discuss  this  in  more 
detail  below  and  in  Section  3.4.5. 

The  value  of  Q/I  for  the  data  of  Jan.  15,  1990  and  Jan.  29,  1990  is  determined 
differently  from  the  way  it  was  determined  for  the  Feb.  4,  1988  data  (for  which  we 
used  Eq.  8).  Since 

Q  _  I\\  ~  Al  _  I\\/h  ~  1 

I  I\\  +  It  I\\/h  +  1 

then  Q/I  can  be  determined  from  the  slope  of  a  plot  of  7X  versus  Jn.  Such  a  plot  is 
shown  for  two  lidar  returns  of  Jan  29,  1990  in  Fig.  10. 

On  this  log-log  plot  the  value  of  I\\/Ix  is  the  vertical  distance  between  the 
data  and  the  line  7j|  =  7X.  For  the  actual  data  analysis  we  determine  the  value 
by  calculating  the  slope  of  the  least  squares  fit  to  the  (non-logged)  data.  The 
uncertainty  is  then  the  standard  deviation  of  the  slope  of  the  least  squares  fit. 

Note  the  fairly  constant  slope  up  to  approximately  1  km.  Beyond  1  km  the 
average  (least  squares)  slope  is  approximately  the  same  but  the  uncertainty  increases 
due  to  greater  uncertainty  in  the  data. 

3.4.3  Range  of  Validity  of  Lidar  Data 

We  expect  the  degree  of  polarization  to  begin  to  decrease  with  range  because  mul¬ 
tiple  snowflake  scattering  should  become  more  dominant.  Although  we  do  observe 
this  decrease  for  many  sets  of  lidar  runs  the  decrease  usually  occurs  beyond  the 
range  that  we  feel  the  lidar  data  is  valid. 

There  is  a  danger  in  deriving  quantities  from  lidar  data  that  are  too  far  out  in 
range.  As  range  increases  signal-to-noise  of  course  decreases.  In  addition,  as  signal 
decreases,  the  discreteness  of  the  data  imposed  by  the  data  acquisition  system  adds 
a  different  type  of  error. 
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Stokes  parameters  for  artificial  data 


Figure  11:  Stokes  parameters  for  four  simulated  lidar  returns. 
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We  have  checked,  by  computer  simulation  of  the  lidar,  the  validity  of  calculat¬ 
ing  degree  of  polarization  from  the  lidar  data  at  far  ranges.  Figure  11  shows  the 
Stokes  parameters  for  four  simulated  lidar  returns  corresponding  to  the  four  po¬ 
larization  measurements.  These  data  were  manufactured  by  making  the  following 
assumptions: 

•  The  data  obeys  the  single  scatter  lidar  equation. 

•  The  extinction  coefficient  is  2.1  km-1,  independent  of  range  (similar  to  typical 
lidar  runs  during  moderate  snow). 

•  The  backscatter  coefficient  is  independent  of  range. 

•  The  degree  of  polarization  is  0.25,  independent  of  range. 

•  The  intensity  of  the  polarization  component  of  the  lidar  return  parallel  to 
that  of  the  laser  is  4000  (transient  recorder  units)  at  its  peak.  The  transient 
recorder  data  word  is  12  bits,  corresponding  to  a  dynamic  range  of  4096.  The 
photomultiplier  gain  function  is  included  in  the  analysis,  which  effectively 
increases  the  dynamic  range  (see  Part  II). 

•  The  intensity  of  the  other  polarization  components  are  chosen  to  be  in  the 
same  proportion  as  one  particular  lidar  run  for  which  the  degree  of  polarization 
is  0.25. 

•  Each  individual  return  from  200  laser  shots  has  Poisson  noise  added.  In 
addition  a  background  level  with  its  own  Poisson  noise  is  added.  The  returns, 
with  the  noise,  are  then  digitized  and  added  to  arrive  at  the  total  accumulated 
return. 

Figure  11  indicates  that  the  degree  of  polarization  is  valid  out  to  the  range  where 
the  signal  minus  the  background  becomes  less  than  1  bit.  This  is  true  for  a  wide 
range  of  degrees  of  polarization,  signal  levels,  and  background  levels,  as  long  as 
enough  shots  are  averaged  over  so  that  the  noise  fluctuations  are  low  enough.  The 
discreteness  imposed  by  digitization  is  the  deciding  factor.  This  is  be  expected  since 
noise  fluctuations  are  merely  fluctuations  about  the  “true”  level.  Taking  enough 
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shots  (if  possible)  should  average  this  effect  out.  However,  if  a  signal  is  less  than  1 
bit,  and  its  fluctuation  is  smaller  than  the  amount  that  would  bring  it  up  to  1  bit, 
then  zero  signal  will  be  recorded,  independent  of  number  of  laser  shots. 

The  degree  of  polarization  and  Stokes  parameters  that  we  have  shown  in  this 
report  corresponds  to  the  portions  of  lidar  returns  which  are  greater  than  1  bit  after 
the  background  is  subtracted  out. 

3.4.4  Snowflake  Scattering  Scenario 

We  present  a  scenario  for  snowflake  scattering  based  on  the  data  and  conclusions  of 
Section  3.4.2.  We  divide  snowflake  scattering  into  two  regimes:  1)  single  snowflake 
scattering  (which  we  observe  and  have  reported  on  above)  and  2)  multiple  snowflake 
scattering  (which  we  do  not  observe  because  of  the  lack  of  validity  of  the  data  at  far 
ranges).  In  the  first  regime  multiple  snowflake  scattering  presumably  occurs  but  is 
negligible.  In  the  second  regime  multiple  snowflake  scattering  is  dominant. 

The  single  snowflake  scattering  regime  can  be  viewed  as  a  multiple  scattering 
phenomenon.  From  the  point  of  view  of  a  photon,  each  individual  snowflake  is  a  col¬ 
lection  of  surfaces  with  random  orientation.  A  photon  enters  the  snowflake  region, 
reflects  from  or  transmits  through  one  or  more  of  the  surfaces,  then  leaves.  Since 
snowflakes  are  large  compared  to  the  wavelength  each  reflection  and  transmission 
should  be  described  fairly  well  by  the  Fresnel  reflection  and  transmission  coefficients 
(discussed  in  more  detail  in  the  theory  Section  3.4.5). 

Each  reflection  or  transmission  alters  the  polarization  of  the  photon  in  a  well 
defined  way.  Some  interactions  do  not  alter  the  polarization  (e.g.,  a  single  back¬ 
ward  reflection  of  linearly  polarized  incident  light).  Because  of  the  large  number 
of  photons  entering  the  snowflake  region,  and  because  of  the  randomness  of  the 
orientations  of  the  surfaces,  we  expect  that  any  photon  changing  its  polarization 
state  becomes  one  of  many  photons  in  one  of  the  many  polarization  states  different 
from  the  incident  photon  polarization  state.  These  photons  contribute  to  the  unpo¬ 
larized  portion  (random  polarization)  of  the  backscattered  radiation.  The  photons 
that  do  not  change  their  polarization  state  contribute  to  the  polarized  portion  of 
the  radiation. 
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Since  this  is  statistical,  we  expect  the  degree  of  polarization  to  have  a  spread 
about  some  value.  The  spread  should  be  smaller  for  greater  numbers  of  snowflakes 
contributing  to  the  backscatter  signal.  In  the  single  snowflake  scatter  regime  the 
received  backscattered  signal  is  the  sum  of  backscattered  signals  from  the  many 
snowflakes  in  a  viewing  volume.  This  may  be  quite  large.  For  example,  for  a 
digitization  time  of  0.1  nsec  and  a  range  of  200  m,  the  receiver  viewing  volume 
is  2.9  m8.  When  the  snow  density  is  0.015  cm-8  (midrange  for  the  plots  shown) 
there  are  43500  snowflakes  in  this  volume.  The  total  number  of  snowflakes  per 
viewing  volume  to  be  considered  is  this  number  times  the  number  of  laser  shots 
being  averaged  over  (200  for  these  data). 

This  number  is  large  enough  that  we  expect  the  degree  of  polarization  to  be  quite 
insensitive  to  snow  density.  However,  when  the  snow  density  is  too  low  scattering 
from  other  particles  that  are  present  besides  snowflakes  becomes  more  prominent. 
If  these  particles  cause  the  radiation  to  preserve  its  polarization  upon  backscatter, 
then  the  degree  of  polarization  increases. 

We  must  add  a  further  comment  due  to  a  characteristic  of  the  plot  in  Fig.  8.  That 
plot  shows  Q//  versus  lidar  run  (for  all  the  runs  of  Feb.  4,  1988)  for  different  ranges. 
Each  range  corresponds  to  a  particular  digitization  by  the  transient  recorders.  The 
plot  indicates  that  the  degree  of  polarization  increases  with  range  up  to  510  m.  (For 
ranges  greater  than  510  m,  which  we  do  not  show,  the  degree  of  polarization,  when 
it  can  be  calculated,  is  fairly  constant  and  approximately  equal  to  the  510  m  value 
for  each  lidar  run.) 

We  do  not  have  a  definite  explanation  for  this  observation.  We  hypothesize 
that  the  single  snowflake  phase  function  for  the  polarized  portion  of  the  radiation 
is  highly  peaked  in  the  180°  direction,  whereas  the  unpolarized  portion  is  fairly 
isotropic.  This  is  a  reasonable  assumption  since  the  180°  direction  is  almost  en¬ 
tirely  composed  of  radiation  that  reflects  off  a  snowflake  surface  that  happens  to 
be  oriented  perpendicular  to  the  direction  of  propagation  of  the  incident  radia¬ 
tion.  Therefore,  as  range  increases  the  receiver  views  a  smaller  portion  of  the  phase 
function  centered  in  the  180°  direction,  or  a  higher  proportion  of  radiation  that  is 
polarized. 
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3.4.5  Theory 

The  theoretical  issues  associated  with  the  scenario  described  above  can  be  separated 
into  three  parts: 

1.  single  scattering  from  individual  snowflakes  (the  fundamental  interaction), 

2.  the  sum  of  many  single  scatters  in  a  receiver  viewing  volume,  and 

3.  multiple  snowflake  scattering. 

Each  part  requires  the  solution  to  the  previous  part,  although  each  can  be  for¬ 
mulated  indpendently.  We  have  addressed  parts  1  and  2  for  snow  scattering.  We 
discuss  these  below.  We  have  only  addressed  part  3,  multiple  scattering,  for  spher¬ 
ical  particles.  We  discuss  multiple  scattering  in  Section  6. 

Part  2:  sum  of  many  single  scatters 

We  discuss  part  2  first  since  we  actually  did  this  work  first.  We  developed  a 
Monte  Carlo  type  model  to  predict  the  degree  of  polarization  of  the  sum  of  radiation 
backscattered  from  the  many  snowflakes  in  a  single  lidar  volume.  The  scattered 
radiation  from  an  individual  snowflake  is  assumed  to  be  in  a  definite,  although 
random,  polarization  state.  Since  we  did  not  have  a  model  for  single  snowflake 
scattering,  we  had  to  assume  a  probability  of  final  states  in  a  certain  way,  which 
we  describe  below.  Adding  the  Stokes  vectors  from  all  the  snowflakes  in  a  viewing 
volume  (see  Eq.  7),  and  making  the  calculation  many  time,  then  gives  a  probability 
for  the  degree  of  polarization  for  the  lidar  return. 

The  calculation  is  performed  in  the  following  way: 

•  For  each  of  N  snowflakes,  choose  ex  and  e„,  the  intensity  of  the  backscattered 
radiation,  from  each  snowflake,  polarized  in  the  x  and  y  directions  respec¬ 
tively.  Choose  by  assuming  each  amplitude  to  be  random  variables  uniformly 
distributed  on  the  interval  [0,1]. 

•  For  each  of  N  snowflakes,  choose  the  phase  difference  A  between  these  two 
amplitudes  by  assuming  it  to  be  a  random  variable  uniformly  distributed  on 
the  interval  {-ir,7rj. 
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•  For  each  of  N  snowflakes,  calculate  the  Stokes  parameters. 

•  Choose  a  number  M  of  scatterers  which  preserve  the  polarization  state  of  the 
incident  radiation  (coherent  scatterers).  The  Stokes  vector  for  the  backscat- 
tered  radiation  from  each  of  these  scatterers  is  [1,1, 0,0]. 

•  Calculate  the  degree  of  polarization  of  light  backscattered  from  snowflake  and 
coherent  scatterers  in  the  following  manner: 

[(m + e"i  «.)’  +  (e~  ,  u.)2  +  (z!i,  v<)‘ 

M  +  EX,  /, 

•  Repeat 

Figure  12  shows  the  result  of  a  Monte  Carlo  calculation,  performed  on  an  80286 
based  computer,  for  10  snowflakes  and  0  coherent  scatterers.  The  figure  shows  plots 
of  the  probability  density  functions  for  the  degrees  of  polarization,  linear  polariza¬ 
tion,  and  circular  polarization.  The  program  was  stopped  when  the  probability 
curves  looked  reasonably  smooth. 

The  degree  of  circular  polarization,  which  is  only  shown  for  postive  values,  is  ac¬ 
tually  symmetric  about  zero.  This  is  because  of  the  assumption  that  the  probability 
of  the  phase  difference  between  ez  and  e„  is  symmetric  about  zero  phase. 

Figures  13a  and  13b  show  two  plots  which  are  a  compendium  of  the  information 
from  many  individual  Monte  Carlo  calculations  of  type  shown  in  Fig.  12.  Figure  13a 
is  a  plot  of  the  peaks  of  the  probability  curves  for  the  degree  of  polarization,  as  a 
function  of  the  ratio  of  coherently  scattered  power  to  power  scattered  by  snowflakes, 
for  different  number  of  snowflakes  in  the  viewing  volume.  The  peak  of  each  proba¬ 
bility  curve  is  approximately  equal  to  the  average  of  the  probability  curve.  The  plot 
in  Fig  13b  is  similar,  except  that  it  shows  the  standard  deviation  of  each  probability 
curve. 

If  we  choose  one  snowflake  and  zero  coherent  scatterers  in  the  volume  then  there 
is  a  probability  of  one  that  the  degree  of  polarization  is  one  (since  each  scatter  al¬ 
ways  results  in  a  definite  polarization  state).  Increasing  the  number  of  coherent 
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Monte  Carlo  results  for 
10  snow  particles 
0  coherent  scatterers 


□  deg  of  pol 

x  deg  of  lin  pol 

v  deg  of  circ  pol 


degree 


Figure  12:  Monte  Carlo  calculation  for  10  snowflakes  and  0  coherenent  scatterers. 
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Coherent  power/Snow  power 


Coherent  power/ Snow  power 


Figure  13:  a)  Average  and  b)  standard  deviation  of  the  degree  of  polarization  versus 
ratio  of  the  number  of  snow  to  coherent  scatterers. 
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scatterers  for  a  single  snowflake  at  first  decreases  the  most  probable  value  of  de¬ 
gree  of  polarization  because  a  coherent  scatter,  by  definition,  always  result  in  one 
particular  polarization  state  (represented  by  the  Stokes  vector  [1,1, 0,0]). 

As  the  number  of  snowflakes  is  increased,  the  most  probable  value  of  degree  of 
polarization  decreases.  However,  the  standard  deviation  decreases  as  well. 

As  the  number  of  snowflakes  approaches  infinity  a  plot  of  the  most  probable  val¬ 
ues  of  degree  of  polarization  as  a  function  of  the  ratio  of  coherent  to  snow  backscat- 
tered  power  approaches  the  curve  labeled  “oo”  on  Fig.  13a.  This  result  is  expected 
due  to  the  form  of  Eq.  6  which  we  rewrite  here: 

S  —  =  ^p/^u 

Ip  + 1 u  Ip/Iu  +  1 

The  most  probable  value  of  degree  of  polarization  for  an  infinite  number  of  snowflakes 
and  zero  coherent  scatterers  is  zero.  That  is,  all  the  backscattered  radiation  from 
snow  is  randomly  polarized.  Therefore,  Ju  is  all  the  radiation  backscattered  from 
snow.  Ip  is  all  the  radiation  backscattered  from  the  coherent  scatterers  since  this 
radiation  is  polarized  and  in  the  same  state.  The  ratio  Iu/Ip  for  an  infinite  number 
of  snowflakes  is  therefore  the  ratio  of  snow  to  coherent  backscattered  power,  the 
variable  on  the  horizontal  axis  of  Fig.  13a.  The  “oo”  curve  in  that  figure  is  simply 
a  plot  of  6  versus  /u//p  from  Eq.  10. 

We  see  that  for  the  many  thousands  of  snowflakes  in  a  lidar  receiver  viewing 
volume  (calculated  in  the  previous  section),  it  is  valid  to  use  the  “oo”  curve. 

Part  1:  single  snowflake  scattering 

The  results  of  a  single  snowflake  scattering  theory  can,  in  principle,  be  used 
to  determine  the  actual  probability  density  functions  for  ex,  cv,  and  A  which  are 
needed  for  the  sum-of-single-scatters  model  given  above.  Recall  that  we  assumed 
these  random  variables  to  be  uniformly  distributed  on  some  range. 

We  have  developed  a  preliminary  model  for  single  snowflake  scattering,  although 
we  did  not  apply  it  to  the  sum-of-single-scatters  model.  As  we  determined  above, 
the  large  number  of  snowflakes  in  a  viewing  volume  allows  us  to  assume  that  the 
degree  of  polarization  lies  on  the  “oo”  curve  shown  in  Fig.  13a.  This  assumption 
is  valid  as  long  as  we  choose  probability  density  functions  such  that  the  average 
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degree  of  polarization  decreases  as  snowflake  density  increases. 

Our  model  for  single  snowflake  scattering  is  based  on  the  scenario  we  presented 
in  Section  3.4.2  on  page  24.  The  following  list  summarizes  the  scenario  and  the 
assumptions  used  by  the  model: 

•  Light  that  hits  a  snowflake  undergoes  a  series  of  one  or  more  transmissions 
and/or  reflections  from  randomly  oriented  surfaces  of  the  snowflake. 

•  The  direction  of  the  normal  vector  of  one  of  these  surfaces  is  a  random  variable 
distributed  uniformly  in  solid  angle  (i.e.,  P(0,  <f>)dCl  =  l/4ndCl). 

•  Each  interaction  (reflection  or  transmission)  is  described  by  Snell’s  law  and 
the  Fresnel  reflection  and  transmission  coefficients. 

•  After  a  given  interaction  all  the  light  which  was  reflected  or  transmitted  is 
available  for  the  next  interaction. 

•  Light  that  changes  polarization  is  considered  part  of  the  unpolarized  portion 
of  the  backscattered  light  (/„).  Light  that  has  polarization  preserved  is  part 
of  the  polarized  portion  (/p).  The  degree  of  polarization  from  the  snowflake 
is  then  6  =  IP/(IU  +  Jp). 

Consider  only  the  set  of  interactions  that  attenuate  light  the  least  (see  Figs.  14a 
through  14c): 

1.  1  reflection  (preserves  polarization), 

2.  2  reflections  (may  change  polarization),  and 

3.  1  transmission-1  reflection-1  transmission  (may  change  polarization). 

An  additional  reflection  will  attenuate  the  light  further  by  an  amount  on  the  order 
of  the  average  Fresnel  reflectivity  of  0.12  (averaged  over  angle  of  incidence  and 
polarization). 

We  sum  up  the  amount  of  polarized  and  unpolarized  power  from  each  of  these 
interactions.  Since  the  degree  of  polarization  is  only  dependent  on  the  ratio,  we  can 
assume  an  incident  power  of  one. 
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Figure  14:  The  three  interactions  considered  for  the  single  snowflake  scatter  model. 
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For  the  first  two  interactions  above  we  analytically  calculated  the  polarized  and 
unpolarized  backscattered  power.  For  the  third  interaction  we  wrote  a  computer 
program  on  an  80386  based  computer.  An  additional  approximation  we  made  in 
regard  to  the  third  interaction  is  discussed  below. 

The  backscattered  powers  from  each  interaction  are  summarized  below,  a  is  the 
receiver  field  of  view. 


interaction 

h 

Iu 

1 

6.1  xl0~4  a2 

0 

2 

4.6x  10-4  a2 

4.2X10"6  a2 

3 

1.6x10-®  a2 

8.0X10'4  a2 

From  this  we  calculate  the  degree  of  polarization: 


Ip  _  6.1  x  10~4  +  4.6  X  10~4  +  1.6  x  10~®  _ 

Tu  ~  4.2  x  lO-6  +  8.0  x  10~4  "* 


6  = 


v/« 


=  0.57 


in) 


(12) 


1  +  Jp/Ju 

The  result  of  0.57  for  the  degree  of  polarization  is  not  far  off  from  the  0.3-0. 5  degree 
of  polarization  that  we  experimentally  observe  for  single  snow  scattering  (see  Fig.  8 
on  page  18).  That  such  an  approximate  theory  would  give  a  result  within  an  order 
of  magnitude  is  an  indication  that  the  basic  assumed  interactions  must  have  some 
validity. 

We  must  qualify  our  result  for  the  intensity  of  the  unpolarized  radiation  for 
the  third  interaction  given  above.  The  results  of  the  computer  calculation  actually 
give  a  value  of  zero  for  Iu.  The  radiation  that  is  scattered  back  to  the  receiver  is 
completely  polarized.  A  large  majority  of  the  radiation  (0.036a2)  is  totally  internally 
reflected  at  the  third  surface,  and  must  go  on  to  either  interact  with  another  surface 
or  leave  the  snowflake  region.  We  assumed  that  this  radiation  is  an  unpolarized, 
isotropic  source  of  radiation  for  one  addition*  ’  reflection.  Summing  up  all  of  this 
radiation  that  can  get  back  to  the  receiver,  making  the  same  assumption  as  before 
about  the  random  orientation  of  the  surfaces,  gives  the  above  result  for  7U  for  the 
third  interaction. 


33 


4  Mie  Scattering 

4.1  Introduction 

Mie  scattering  is  the  overall  guiding  theoretical  mechanism  for  scattering  in  our 
lidar  work.  Although  it  is  not  always  applicable,  it  is  usually  the  best  theory  which 
is  exactly  soluble. 

Early  on  in  this  program  we  studied  Mie  scattering,  from  the  point  of  view  of 
understanding  its  validity  in  our  lidar  experiments.  In  particular,  we  wanted  to 
determine  the  validity  of  the  widely  used  assumption  that  the  extinction  coefficient 
is  proportional  to  the  backscatter  coefficient  for  a  collection  of  Mie  scatterers  (e.g., 
aerosols)  with  a  size  distribution  typical  of  the  atmosphere. 

We  developed  computer  programs  which  calculate,  for  a  range  of  particle  sizes, 
the  total  Mie  scattering  cross  section  and  the  differential  scattering  cross  section 
for  a  180°  scattering  angle,  and  convolve  these  cross  sections  with  arbitrarily  chosen 
size  distributions.  Results  are  presented  in  Section  4.4.  Section  4.2  is  a  background 
discussion  of  Mie  scattering.  It  presents  some  theoretical  background  and  concludes 
with  expressions  for  the  Mie  scattering  total  cross  section  and  differential  scattering 
cross  section  for  180°.  Section  4.3  describes  our  computer  model. 

4.2  Theoretical  Background 

Mie  scattering  refers  to  the  scattering  of  a  monochromatic  electromagnetic  plane 
wave  by  a  spherically  shaped,  homogeneous,  isotropic  dielectric  and  conducting 
medium  embedded  in  an  infinite,  homogeneous,  isotropic,  dielectric,  nonconduc- 
tiong  medium.  The  analytic  solution  of  this  problem  was  published  by  Gustav  Mie 
in  1908  [3].  It  is  the  solution  of  Maxwell’s  equations  with  boundary  conditions 
corresponding  to  the  above  scenario. 

The  electric  field  of  the  incident  plane  wave  is  assumed  to  be  exp  («(fcx«  —  u >*)]£, 
where  k\  =  k0n j,  is  the  real- valued  index  of  refraction  of  the  infinite  medium, 
and  k0  =  w/c  is  the  vacuum  wave  vector  amplitude.  The  amplitude  of  the  incident 
field  is  unity.  Its  polarization  is  in  the  x  direction.  The  index  of  refraction  nj  of 
the  sphere  is  complex. 

The  scattered  electromagnetic  field  is  the  sum  of  two  subfields.  One  subfield  has 
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a  magnetic  field  with  zero  radial  component  and  the  other  subfield  has  an  electric 
field  with  zero  radial  component.  Each  subfield  is  an  infinite  series  solution  to  the 
equations  and  boundary  conditions.  The  solution  is  cast  in  spherical  coordinates. 

The  electric  and  magnetic  fields  of  the  first  subfield  (zero  radial  magnetic  field), 
designated  with  a  superscript  *,  are  [2] 

1  cos  <j> 


Eer 


El  =  - 


£/(/  +  l)B,eft(fcir)Pu(cos0) 
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The  electric  and  magnetic  fields  of  the  second  subfield  (zero  radial  electric  field), 
designated  with  a  superscript  m  are  [2] 
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In  Eqs.  13  and  14 


is  an  associated  Legendre  polynomial  with  m  =  1, 

Si(p)  =  p\ji{p)  -  »yi(p)], 

ji  and  yi  are  the  usual  spherical  Bessel  functions, 

“prime”  (')  denotes  differentiation  with  respect  to  argument,  and 
the  time  dependence  of  each  field  component  is  e~,ut. 
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The  coefficients  B‘  and  B/*,  which  are  only  functions  of  the  permitivities  and 
conductivity  of  the  sphere  and  surrounding  medium,  are  given  by 

D(  -i+i  (2i  +  1)  \nxl>\{q)Mn<l)  ~  MM{m)  1 


where 


/(/  +  1)  [  nf/(g)0i(ng)  -  {i[q)tpl[hq) 
[21  +  1)  [n^i(g)0|(ng)  -  4>\{q)il>i[hq) 


/(/  +  !)  [  hfi(q)t/>i( 


j(nq)  -  t/>!(q)r/>i(nq)' 
l(nq)  ~  C/(q)^i(nq)  . 


„  n2  e2  Attk 

n  =  —  = - 1- 1 - 

n j  If!  u>ci . 


q  =  kid  (“size  parameter”), 

a  is  the  sphere  radius, 

ci  is  the  permitivity  of  the  infinite  medium, 

c2  is  the  real  part  of  the  permitivity  of  the  sphere, 

k  is  the  conductivity  of  the  sphere,  and 

ipi(p)  =  pji[p )  (Ricatti-Bessel  function). 

The  differential  scattering  cross  section  is  defined  as 

da  I  _  dP,cat[6)/dn 


where 
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(time  averaged  power  scattered  per  unit  solid  angle  into  direction  0 ),  and 

Sine  =  ni^-z  (19) 

57T 

(time  averaged  incident  flux). 

The  total  scattering  cross  section  is  the  integral  of  Eq.  18  over  all  solid  angles. 
This  is  a  very  difficult  task  to  perform,  given  the  form  of  the  scattered  fields.  It 
turns  out  that,  for  large  r  and  for  polarized  incident  wave,  the  total  cross  section 
(scattering  plus  absorption)  can  be  determined  from  the  differential  scattering  cross 


section  in  the  0  =  0  direction  [2j.  That  this  is  the  case  is  due  to  conservation  of 
energy.  The  total  amount  of  radiation  lost  to  scattering  and  absorption  is  just 
the  amount  of  radiation  left  over  propagating  in  the  direction  of  the  incident  wave 

(0  =  o). 

The  total  cross  section,  in  the  context  of  the  current  problem,  is 

o  =  oteat  +  oai,  =  ^Im  ji  •  E{z) ]  (20) 

—4 

where  E(z)  is  the  scattered  electric  field  (E(z))  with  the  term  e'*‘r/r  factored  out, 
and  evaluated  in  the  direction  of  the  propagation  of  the  incident  wave  (z).  The 
dot  product  is  then  taken  with  the  unit  vector  in  the  direction  of  the  incident 
polarization  x.  We  note  that  E  is  easily  extracted  from  Eqs.  13  and  14  because 
fj(fcir)  and  f/(fcjr)  both  go  like  e,k‘r  for  large  k\r  [2).  Equation  20  is  known  as  the 
optical  theorem. 

By  inserting  Eqs.  13  and  14  into  Eq.  20  the  total  cross  section  becomes  [2] 

1  °° 

*  =  «Re  £  H),+‘  n‘  +  i)  (b;  +  bd  .  (21) 

"i  Lt=i 

We  also  write  here  the  differential  scattering  cross  section  for  180°,  which  we 
refer  to  as  the  backscatter  cross  section.  We  will  use  this  in  the  next  section.  The 
backscatter  cross  section  is  determined  by  inserting  Eqs.  13  and  14  into  Eq.  18  and 
setting  0  =  n: 


where  Pt  is  the  Ith  Legendre  polynomial. 

4.3  Computer  Model 

We  have  developed  a  computer  program  which  calculates  the  total  scattering  cross 
section  (Eq.  21)  and  the  backscatter  cross  section  (Eq.  22)  for  a  dielectric,  noncon- 
ductive  sphere  of  arbitrary  size.  For  zero  conductivity  ( k  =  0)  there  is  no  absorption 
and  the  total  cross  section  is  equal  to  the  scattering  cross  section. 

We  use  the  results  of  this  program  to  calculate  the  extinction  coefficient  and  the 
backscatter  coefficient  for  a  collection  of  dielectric,  nonconductive  spheres.  When 
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the  spheres  are  randomly  distributed  and  the  spacing  between  spheres  is  large  com¬ 
pared  to  the  wavelength  the  extinction  coefficient  a  and  b&ckscatter  coefficient  (3 
are 

{;}-f  "“>{(*).}*  {23) 

where  N(a ),  the  size  distribution,  is  the  number  of  spheres,  per  unit  volume,  with 
radii  between  a  and  a  +  da.  The  normalization  is  /“  N(a)da  =  Nt ,  where  Nt  is  the 
total  density  of  particles.  The  size  distribution  can  be  chosen  arbitrarily. 

All  programs  and  subroutines  were  developed  and  run  on  an  80286  based  com¬ 
puter.  The  infinite  sums  in  Eqs.  21  and  22  were  truncated  at  the  first  term  which 
was  less  than  some  fraction  of  the  accumulated  sum  of  the  previous  terms.  This  was 
adequate  for  calculating  cross  sections  for  particles  with  radii  up  to  approximately 
50  times  the  wavelength. 

4.4  Results  and  Conclusions 

We  first  show  results  of  calculations  of  total  and  backscatter  cross  sections  for  indi¬ 
vidual  spheres.  These  results  serve  to  verify  that  the  program  is  working  correctly. 
We  then  show  results  of  calculations  of  extinction  and  backscatter  coefficients  for 
collections  of  spheres  with  size  distributions  typical  of  atmospheric  aerosols. 

4.4.1  Individual  Spheres 

Figure  15  shows  the  total  and  backscatter  cross  section  efficiencies  for  a  single 
sphere  versus  size  parameter.  The  cross  section  efficiency  (total  or  backscatter)  is 
the  cross  section  divided  by  the  cross  sectional  area  7ra2  of  the  sphere,  where  a  is 
the  particle  radius.  The  size  parameter  is  2na/Xi  =  fcjo.  The  indices  of  refraction 
for  these  plots  are  nj  =  1  (surrounding  medium:  vacuum)  and  nj  =  1.33  (sphere: 
water). 

We  note  several  observations  from  these  plots  which  serve  to  verify  the  proper 
operation  of  the  program: 

•  For  small  size  paramater  both  cross  section  efficiencies  go  as  the  particle  size 
parameter  raised  to  the  fourth  power.  For  small  size  parameter  we  expect  the 
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Mie  Scattering  cross  section  efficiencies 
n\  =  1.0  n  2  =  1.33 


Figure  15:  Mie  total  and  backscatter  cross  section  efficiencies  versus  size  parameter 
for  ni  =  1,  n2  =  1.33,  and  k  =  0. 
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Mie  scattering  solution  to  reduce  to  Rayleigh  scattering  solution.  In  Rayleigh 
scattering  the  cross  section  goes  as  (particle  volume)2 /A4,  or  as  rfi/A4.  There¬ 
fore,  cross  section  efficiency  goes  as  r4/ A4,  which  we  observe. 

•  The  ratio  of  backscatter  cross  section  efficiency  to  total  cross  section  efficiency 
to  is  0.12  for  small  size  parameters.  For  Rayleigh  scattering  this  ratio  is 
3/8 ir  «  0.12. 

•  The  total  cross  section  efficiency  approaches  the  value  2  for  large  size  pa¬ 
rameters.  For  large  size  parameters,  we  expect  the  Mie  scattering  solution 
to  reduce  to  the  result  obtained  from  scalar  diffraction  theory.  For  an  arbi¬ 
trary  shaped  obstacle  which  is  large  compared  to  the  wavelength  and  does 
not  transmit  much  light  the  main  contribution  to  the  forward  scattered  light 
is  from  Fraunhofer  diffraction.  For  distances  far  from  the  obstacle  this  total 
cross  section  is  2D,  where  D  is  the  cross  sectional  area  of  the  obstacle  as 
“seen”  by  the  incident  plane  wave.  The  efficiency  is  therefore  2. 

Note  the  very  fast  fluctuations  in  both  cross  sections  as  a  function  of  size  pa¬ 
rameter  for  size  parameters  greater  than  approximately  one.  (Because  of  the  log 
scale  the  fast  fluctuations  in  total  cross  section  efficiency  appear  small.)  This  is  the 
result  of  adding  coherently  all  the  partial  waves  of  each  subfield  (each  term  of  the 
sums  in  Eqs.  21  and  22).  The  ratio  of  backscatter  cross  section  to  total  cross  sec¬ 
tion  in  this  regime  is  very  complicated  and  would  not  be  very  useful  experimentally. 
Fortunately,  a  size  distribution  with  a  variance  large  compared  to  these  fluctuations 
tends  to  average  out  the  fluctuations.  This,  as  we  show  below,  is  what  happens  for 
atmospheric  aerosol  scattering. 

4.4.2  Collection  of  Spheres 

There  are  many  different  size  distributions  which  describe  the  many  different  types 
of  particles  of  the  lower  troposphere.  We  have  not  performed  a  study  of  particle 
sizes  in  the  atmosphere.  Our  intention  here  is  to  consider  size  distributions  generally 
accepted  as  describing  some  atmospheric  constituents  and  to  use  them  with  the  Mie 
formulation  to  determine  the  behavior  of  the  ratio  of  backscatter  coefficient  to  total 
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extinction  coefficient.  We  present  two  size  distributions  below. 

A  size  distribution  used  to  describe  the  fine  water  aerosols  which  form  so-called 
continental  haze  [4]  is 

<»> 

where  i/p  is  a  constant  which  is  usually  between  2.0  and  5.0  and  J  is  a  normalization 
constant.  This  power  law  model  was  originally  propose  by  Junge  [5]. 

A  size  distribtion  that  is  found  to  represent  clouds  and  fogs  [4]  takes  the  form 

AT(fl)  =  C  (— VVa/a-  (25) 

VOm/ 

where  am,  the  mode  radius,  is  the  value  of  a  for  which  N(a)  is  maximum  and  C  is 
the  normalization  constant.  i/p  determines  the  spread  of  the  distribution  about  the 
maximum.  This  model  has  been  discussed  in  detail  by  Deirmendjian  [6,7]. 

Figure  16  shows  a  plot  of  backscatter  efficiency  versus  extinction  efficiency  for 
the  haze  model.  Each  point  represents  a  calculation  of  Equation  23  with  a  different 
value  of  i/p  in  the  continental  haze  size  distribution.  The  values  of  up  are  in  the 
range  of  2-5  in  steps  of  0.05.  The  single  particle  Mie  scattering  formulation  with 
conditions  identical  to  those  of  Fig.  15  are  used. 

Figure  17  shows  four  plots  of  backscatter  efficiency  versus  extinction  efficiency 
for  the  cloud/fog  model.  Each  plot  corresponds  to  a  different  value  of  i/p,  as  marked. 
Each  point  in  each  plot  represents  a  calculation  of  Equation  23  with  a  different  value 
of  kiam  in  the  fog/cloud  size  distribution.  The  values  of  kiam  are  in  the  range  of 
0.1-50.  The  step  is  0.1  in  the  range  0.1-1  and  1  in  the  range  1-50.  Again,  the  single 
particle  Mie  scattering  formulation  with  conditions  identical  to  those  of  Fig.  15  are 
used. 

The  plot  in  Fig.  16  indicates  that  the  assumption  of  constant  backscatter  to 
extinction  ratio  for  the  continental  haze  model  is  fairly  good.  The  plots  of  Fig.  17 
indicate  that  the  assumption  for  the  cloud/fog  model  is  fairly  good  for  mode  radius 
larger  than  a  certain  amount,  dependent  of  vp.  For  large  mode  radius  the  integral  in 
Eq.  23  has  largest  contributions  from  the  scalar  diffraction  regime  of  Mie  scattering. 
For  small  mode  radius  the  assumption  is  not  as  good,  but  still  not  too  bad.  Here 
the  integral  has  large  contributions  from  the  intermediate  regime  between  Rayleigh 
scattering  and  scalar  diffraction. 


41 


backscatter  efficiency 


Figure  16:  Backscatter  versus  extinction  efficiencies  for  different  values  of  i/p  of  the 
continental  haze  model. 
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ter  efficiency  bockscatter  efficiency 


Figure  17:  Backscatter  versus  extinction  coefficients  for  different  values  of  vp  and 
mode  radius  am  of  the  cloud/fog  model. 


5  Extinction 


5.1  Introduction 

A  major  use  of  the  low  altitude  lidar  system  is  in  determining  the  extinction  coef¬ 
ficient  for  a  variety  of  meterological  conditions  in  the  troposphere.  The  extinction 
coefficient  is  determined  from  lidar  data  only  with  the  aid  of  assumptions,  mod¬ 
els,  or  additional  experiments  which  are  independent  of  the  lidar  data.  Our  work 
has  involved  testing  the  validity  of  the  extinction  calculations,  understanding  their 
limitations,  and  enhancing  the  extinction  algorithm. 

Section  5.2  describes  the  standard  Klett  inversion  method  that  is  the  basis  of  all 
of  our  extinction  calculations.  It  describes  the  inversion  algorithm  (Section5.2.1), 
our  methods  of  determining  the  extinction  at  the  final  range  of  interest  (Sec¬ 
tion  5.2.2),  and  some  calculations  with  artificial  data  (Section  5.2.3).  Section  5.3 
shows  the  results  of  extinction  calculations  made  in  clear  air  and  clouds  (Sec¬ 
tion  5.3.1)  and  falling  snow  (Section  5.3.2).  The  results  for  falling  snow  were  used 
for  the  polarization  studies  described  in  Section  3.4. 

5.2  Theoretical  Background 

5.2.1  Inversion 

The  range  resolved  extinction  coefficient  a(r)  is  determined  from  the  lidar  data 
with  the  Klett  inversion  method.  The  Klett  inversion  [8]  is  a  (numerically)  stable 
solution  to  the  single  scatter  lidar  equation.  It  requires  knowledge  of  the  extinction 
coefficient  at  the  final  range  value  of  interest. 

The  single  scatter  lidar  equation  is 

P M  =  P.f  exp  [-2  J'  a(rVr']  (26) 

where  P(r)  is  the  instantaneous  range  resolved  received  power,  P0  is  the  transmitted 
power,  c  the  velocity  of  light,  r  the  pulse  duration,  and  A  the  effective  receiver  area. 
0(r)  is  the  backscatter  coefficient  (km-1sr-1)  and  a(r)  is  the  extinction  coefficient 
(km-1),  both  at  a  range  r  from  the  receiver/transmitter. 

Equation  26  has  two  unknowns:  a(r)  and  /?(r).  To  solve  Eq.  26  we  assume  the 
relationship  (3  =  Cak,  where  C  and  k  are  constants.  We  have  addressed  the  validity 
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(27) 


of  this  relationship  in  Section  4.4.  The  lidar  equation  is  then  rewritten  as 

dS  k  doc  n 

dr  a  dr 


where  5  =  ln[rlP(r)j  is  called  the  log  squared  corrected  power. 

Equation  27  is  a  form  of  Ricatti  equation  which  is  exactly  soluble.  Since  it  is  a 
first  order  differential  equation,  the  solution  is  specified  with  one  arbitrary  initial 
condition.  It  is  typical  to  choose  the  initial  condition  as  the  value  of  the  independent 
variable  (a  in  this  case)  at  the  beginning  of  the  range  interval  for  which  solution  is 
desired.  However,  doing  so  in  this  situation  leads  to  a  numerically  unstable  solution. 

The  numerically  stable  solution  put  forth  by  Klett  involves  choosing  the  “initial” 
condition  at  the  end  of  the  range  interval.  The  solution  is 


exp| 

I 

£  +  f£'«Pl 

F= 

\  dr1 

where  Sf  =  5(r/),  af  =  a[rf),  and  17  is  the  final  range  of  interest.  Equation  28 
is  used  for  r  <  17.  Equation  28  is  the  solution  we  use  for  determining  extinction 
coefficients  from  lidar  data. 

The  Klett  inversion  of  lidar  data  can  be  performed  on-line  with  the  the  lidar 
control  program.  The  details  of  the  implementation  are  described  in  Part  II.  The 
program  will  allow  a  value  for  07  to  be  specified,  or  it  will  choose  on  the  basis  of 
the  one  of  the  algorithms  described  in  the  next  section. 


5.2.2  Initial  Condition 

From  a  mathematical  point  of  view  the  initial  condition  07  is  arbitrary.  There¬ 
fore,  the  only  correct  way  of  specifying  a  is  by  measuring  it.  However,  there  are 
ways  of  making  “smart”  guesses.  Below  we  describe  two  methods  which  we  have 
implemented  in  the  lidar  control  program.  The  first  method,  which  assumes  a  to 
be  constant  in  a  small  interval  near  *7,  is  standard  and  used  throughout  the  lidar 
community.  The  second  method,  which  assumes  a  to  be  a  linear  function  of  r  near 
r/,  is  an  enhancement  that  we  developed. 

Constant  07 
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This  method  assumes  that  a  is  a  constant  function  of  range  near  r/.  The  lidar 
equation  (Eq.  27)  for  constant  a  has  the  solution 

(») 


Therefore,  if  a  is  constant  then  the  slope  of  S  must  be  constant.  The  “best”  a /  can 
therefore  be  found  by  determining  the  linear  least  squares  slope  of  5  for  the  some 
number  of  data  points  near  the  final  range  for  which  a  solution  is  desired. 

For  a  positive  slope  of  S  this  method  breaks  down  because  it  would  result  in  a 
negative  extinction  coefficient.  The  reason  is  that  the  assumption  of  constant  a  is 
no  longer  valid.  This  might  happen  when  the  extinction  coefficient  increases  very 
quickly,  as  is  the  case  at  the  interface  between  clear  air  and  a  cloud.  The  second 
method  of  choosing  a/  does  not  suffer  from  this  problem. 

Linear  a/ 

This  method  assumes  that  a  is  an  arbitrary  linear  function  of  range  (not  neces¬ 
sarily  constant)  near  r/.  If  we  write 

a  =  Ar  +  B,  (30) 

where  A  and  B  are  unknown  constant  parameters  to  be  determined,  then  the  lidar 
equation  (Eq.  27)  becomes 

f  =  I7T5-2^  +  fl>-  <31> 

The  “best”  A  and  B  are  found  by  least  square  fitting  the  data  near  tj  to  Eq.  31. 

That  is,  the  best  A  and  B  are  those  that  minimize  the  function 

=  +  (32) 

where  r<  is  the  range  corresponding  to  the  ith  data  point,  5,  =  <S'(r,),  rn  =  rf,  and 
m  is  one  less  than  the  number  of  data  points  considered  near  rj  for  the  least  squares 
fit. 

Minima  of  the  function  I(A,  B)  with  respect  to  A  and  B  would  typically  be 
found  by  solving 

"  =  0 
dA 

~  =  0.  (33) 
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However,  this  nonlinear  system  of  equations  cannot  be  solved  analytically  and  we 
therefore  resort  to  a  nun-  erical  technique  to  minimize  I. 

A  simplex  method  [9]  is  used  to  minimize  /  with  respect  to  A  and  B.  Initially, 
three  points  are  chosen  in  the  A  —  B  plane.  One  of  these  points  cam  be  A  =  0  and 
B  equal  to  the  a/  found  from  the  first  method  of  determining  a/  (“constant  a/* 
method,  described  above).  The  other  two  points  are  near  the  first  point.  Next, 
of  these  three  points  the  point  which  gives  maximum  I  is  reflected  across  the  line 
connecting  the  other  two  points.  The  value  of  I  at  this  new  point  will  be  less  than 
the  values  at  the  other  points  (for  well-behaved  I).  This  process  is  continued  until 
the  percent  change  of  the  value  of  I  at  the  new  point,  relative  to  the  previous  point, 
is  less  than  some  desired  tolerance. 

5.2.3  Extinction  Calculations  with  Artificial  Data 

We  have  tested  the  extinction  algorithm  with  artificially  created  data.  Assuming 
an  extinction  profile  shown  by  the  solid  line  in  Fig.  18b  we  calculate  the  lidar  power 
which  would  give  this  profile  (Fig.  18a.  We  then  calculate  the  extinction  coefficient 
with  the  Klett  algorithm.  The  data  shown  in  Fig.  18a  could  be  an  idealization  of  a 
lidar  return  from  a  cloud. 

The  nonsolid  lines  in  Fig.  18  show  the  calculated  extinction  for  different  choices 
of  ctf.  The  final  range  is  chosen  past  the  cloud.  Two  of  the  a/  values  are 
chosen  by  the  two  methods  described  in  the  previous  section  (constant  and  linear 
a.j  methods.  The  calculated  extinctions  for  these  two  values  of  aj  lie  almost  entirely 
on  the  model  extinction  curve  and  are  difficult  to  see.  Both  values  of  a/  for  the  two 
methods  are  almost  the  same. 

The  other  values  of  a/  are  chosen  to  illustrate  the  stability  of  the  inversion 
algorithm.  Note  that  values  of  o/  quite  different  from  the  true  a/  give  solutions 
that  quickly  converge  to  the  true  solution. 

Figures  19  and  20  are  similar  to  Fig.  18  except  that  random  Gaussian  noise 
is  added  to  the  artificially  created  lidar  return  before  calculating  the  extinction 
coefficient.  The  standard  deviation  of  the  Gaussian  noise  of  Figure  20  is  ten  times 
that  of  Fig.  19.  In  both  figures  the  extinction  calculation  is  initiated  in  the  cloud, 
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Figure  18:  a)  Artificial  lidar  return,  b)  Solid  line:  model  extinction  profile.  Others: 
Klett  inversions  with  different  a/’s. 
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since  the  data  after  the  cloud  is  too  noisy.  Both  figures  show  calculations  with  the 
linear  and  constant  a/  methods.  The  linear  aj  calculation  is  initiated  in  a  region 
where  the  constant  ctf  method  could  not  be  used  because  dS/dr  >  0. 

5.3  Extinction  Calculations  with  Real  Data 

5.3.1  Clear  Air  and  Clouds 

Figure  21  shows  a  lidar  return  from  clear  air  along  with  a  Klett  inversion.  A  value 
of  a/  was  chosen  by  assuming  a  to  be  constant  over  the  final  ten  range  points 
considered. 

Figure  22  shows  a  lidar  return  with  a  cloud,  and  a  Klett  inversion  through  the 
cloud.  A  value  of  a/  was  chosen  by  assuming  a  to  be  constant  over  the  final  four 
range  points  considered. 

The  lidar  returns  in  these  two  figures  were  acquired  during  daytime  in  Sudbury, 
MA.  We  find  that  the  extinction  coefficient  during  clear  air  conditions  there  to  be 
in  the  range  of  0.1  to  0.5  km-1.  The  extinction  coefficient  in  the  cloud  of  Fig.  22 
indicates  that  the  laser  power  was  too  low  to  observe  the  top  side  of  the  cloud. 

5.3.2  Falling  Snow 

We  find  that  the  the  technique  of  calculating  the  extinction  coefficient  from  lidar  is 
quite  robust.  This  is  revealed  by  our  calculations  using  data  from  falling  snow.  In 
our  polarization  studies  of  Section  3.4  we  presented  and  used  the  lidar  calculated  ex¬ 
tinction  coefficients  in  several  plots.  Here  we  show  the  validity  of  this  by  comparing 
these  extinction  calculations  with  data  from  other  groundbased  instruments. 

During  acquisition  of  lidar  data  in  falling  snow  on  Feb.  4,  1988  data  was  simul¬ 
taneously  acquired  with  a  Particle  Measurement  Systems  (PMS)  rain  distrometer 
and  a  visibility  meter.  Both  instruments  were  located  on  the  ground  approximately 
four  meters  from  the  lidar  transmitter/receiver.  The  visibility  meter  measures  the 
extinction  coefficient  directly.  It  acquires  data  at  a  ten  second  rate. 

The  rain  distrometer  measures  the  snow  fall  rate  (cm-2sec-1)  and  the  distri¬ 
bution  of  snowflake  diameters  in  one  dimension  (snowflakes  interrupt  a  laser  beam 
and  cast  a  shadow  on  a  64  pixel  line  detector).  We  assume  that  each  snowflake  is 
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Figure  19:  a)  Artificial  lidar  return  with  Gaussian  noise,  b)  Dotted  line:  model 
extinction.  Solid  lines:  calculated  extinctions. 
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Figure  20:  a)  Artificial  lidar  return  with  Gaussian  noise,  b)  Dotted  line:  model 
extinction.  Solid  lines:  calculated  extinctions. 
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circular  with  the  measured  diameter.  Then,  the  density  of  particles  Ni  with  radii 
a,  (t  =  1, . . .  ,64;  a*  are  the  measured  diameters  divided  by  two)  is  determined  by 
multiplying  by  the  snow  fall  velocity,  which  we  asuume  to  be  50  cm/sec  [10]  for  all 
snowflakes.  The  total  density  Nt  is  then 

ty^A^Aa,  (34) 

«=i 

where  Aa  is  the  distance  between  pixels. 

The  extinction  coefficient  a  is  determined  by  assuming  an  extinction  efficiency 
of  2  (which  is  the  case  for  particles  large  compared  to  a  wavelength — see  Section  4). 
Then, 

64 

a  =  2  ^2  na^NiAa.  (35) 

«= l 

Figure  23  shows  the  extinction  coefficient  determined  with  the  visibility  meter 
versus  lidar  calculated  extinction  coefficient.  Figure  24  shows  extinction  coefficient 
determined  with  the  rain  distrometer  versus  lidar  calculated  extinction.  Each  point 
in  each  graph  corresponds  to  a  particular  lidar  return.  The  lidar  extinction  shown 
is  the  average  extinction,  calculated  with  the  Klett  inversion,  over  the  interval  180- 
500  m.  Visibility  and  rain  distrometer  data  for  each  point  corresponds  to  an  average 
over  a  four  minute  interval  starting  at  the  time  of  the  lidar  return. 

Also  shown  in  Figs.  23  and  24  are  the  correlation  coefficients  and  the  least  squares 
slopes  between  the  visibility  meter  and  lidar  and  between  the  rain  distrometer  and 
lidar  respectively.  The  correlations  are  0.95  and  0.82  in  the  first  and  second  plots 
respectively.  These  are  fairly  high  correlations  and  we  feel  that  this  indicates  that 
the  lidar  is  a  useful  tool  for  determining  extinction  in  snow. 

There  are  seveal  reasons  why  the  correlations  are  not  unity.  The  lidar  measures 
a  spatially  averaged  extinction  coefficient,  averaged  over  a  relatively  short  time 
interval.  The  visibility  meter  and  rain  distrometer  measure  a  temporally  averaged 
extinction  coefficient,  averaged  over  a  relatively  small  region  of  space.  Furthermore, 
the  regions  of  space  do  not  overlap. 

The  lidar  extinction  calculation  uses  the  assumption  of  constant  ratio  between 
backscatter  and  extinction  coefficients  along  the  entire  propagation  path.  We  expect 
this  to  be  true  for  snow  only  on  average.  Because  of  the  random  shape  and  orien- 
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Figure  23:  Visibility  meter  extinction  coefficient  versus  lidar  extinction  coefficient. 
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Figure  24:  Rain  distrometer  extinction  coefficient  versus  lidar  extinction  coefficient. 
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tation  of  snow,  this  ratio  might  fluctuate  about  some  average  value,  much  like  the 
degree  of  polarization  of  the  backscattered  radiation,  as  we  discussed  in  Section  3.4. 

6  Multiple  Scattering 

6.1  Introduction 

PhotoMetrics  has  done  some  preliminary  modeling  of  multiple  scattering  effects 
for  lidar.  Although  we  do  not  have  results  from  this  work  yet,  we  have  defined  a 
direction  along  which  these  modeling  efforts  could  be  pursued.  We  have  determined 
a  theoretical  approach  and  have  derived  expressions  for  the  backscattered  power 
based  on  this  appoach. 

Section  6.2  briefly  discusses  the  theoretical  background.  It  identifies  the  different 
theoretical  approaches  which  are  used  in  multiple  scattering  theory. 

Section  6.3  discusses  the  radiative  transfer  equation,  the  approach  that  we  have 
chosen  to  take. 

Section  6.4  describes  our  application  of  the  radiative  transfer  equation  to  a 
particular  lidar  problem. 

6.2  Background 

Multiple  scattering  increases  in  significance  as  optical  path  length  increases.  For 
lidar  it  is  important  to  consider  multiple  scattering  after  a  few  hundred  meters  in 
cumulus  clouds,  a  few  kilometers  in  moderate  snow,  and  greater  than  approximately 
ten  kilometers  in  clear  urban  aerosols. 

Many  techniques  exist  which  include  multiple  scattering  effects  in  calculating 
light  propagation  in  a  random  distribution  of  scatterers  such  as  the  atmosphere. 
These  techniques  are  either  analytical  or  Monte  Carlo.  The  analytical  techniques 
are  either  transport  theory  or  multiple  scattering  theory.4  Monte  Carlo  techniques 
are  applied  to  either  of  these  analytical  techniques. 

Transport  theory,  which  treats  the  medium  as  a  continuum,  uses  the  radiative 
transfer  equation  for  the  specific  intensity.  The  radiative  transfer  equation  is  similar 
to  the  Boltzmann  equation  for  the  phase  space  distribution  function  of  particles  in 
4Thi»  terminology  vs  from  Ishimaru  [11]. 
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a  gas  or  plasma.  The  Boltzmann  equation  for  gases  is  typically  a  seven  dimensional 
problem:  three  spatial  dimensions,  three  momentum  or  velocity  dimensions,  and 
time.  The  radiative  transfer  equation  is  a  six  dimensional  problem  for  each  of 
the  four  polarization  degrees  of  freedom:  three  spatial  dimensions,  two  direction 
dimensions,  and  time.6  Thus,  the  radiative  transfer  equation  may  have  as  many  as 
24  degrees  of  freedom. 

Multiple  scattering  theory  treats  the  medium  as  a  collection  of  discrete  scatter- 
ers.  It  starts  with  Maxwell’s  equations,  or  a  reduced  form  of  Maxwell’s  equations 
such  as  the  scalar  wave  equation,  obtains  a  solution  for  single  particle  scattering, 
and  then  introduces  the  interaction  for  many  particles.  This  is  typically  done  by 
summing  over  all  possible  single  scatters  that  can  occur.  Solutions  are  typically  of 
the  form  of  statistical  averages  and  correlations  of  field  quantities. 

We  have  decided  to  formulate  the  multiple  scattering  problem  using  the  radiative 
transfer  equation. 

6.3  Radiative  Transfer  Equation 

The  time  dependent  radiative  transfer  equation  is  four  coupled  integro-differential 
equations  in  the  four  Stokes  parameters.  We  will  only  consider  the  radiative  transfer 
equation  for  the  specific  intensity,  the  first  Stokes  parameter.  Therefore,  we  are 
neglecting  changes  in  the  polarization  of  the  radiation  due  to  multiple  scattering. 

The  time  dependent  radiative  transfer  equation  is  then 

+  «  •  Vr/(r,  s,  t )  +  a</(r,  s,t)  =  ^  f  p(s,  S')/(r,  s',  t)dCl'  (36) 
Ot  47T  J 4r 

where  J(r,  s,t)  is  the  specific  intensity  (W/m2-sr-sec-Hz)  at  position  r  propagating 
in  direction  S  at  time  t,  at  is  the  total  extinction  cross  section  (total  cross  section 
times  particle  density),  a,  is  the  extinction  cross  section  due  to  scattering  (scatter¬ 
ing  cross  section  times  particle  density),  and  p(s,i')  is  the  scattering  phase  function 
from  direction  s  to  s'.  The  phase  function  is  normalized  such  that 

-W  p(s,s')dn'  =  ^-  =  W0.  (37) 

4?r  J ix  Of 

6  The  direction  in  the  radiative  transfer  equation  is  equivalent  to  velocity  in  the  Boltsmann 
equation.  However,  two  directions  completely  specify  the  third  direction  since  the  magnitude  is  the 
speed  of  light. 
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W0  is  called  the  single  particle  albedo.  Equation  36  assumes  no  sources  of  radia¬ 
tion  other  than  those  specified  by  boundary  conditions  (e.g.,  there  is  no  blackbody 
radiation  emitted  at  every  point  in  space). 

For  lidar  it  is  usually  possible  to  reduce  the  number  of  degrees  of  freedom  of 
Eq.  36  to  achieve  a  tractable  problem.  For  example,  for  spherical  particles  and 
cylindrically  symmetric  incident  laser  radiation,  the  specific  intensity  is  independent 
of  azimuth  angle.  For  isotropic  scattering  the  specific  intensity  is  independent  of  the 
angular  direction  orthogonal  to  the  laser  propagation  direction.  For  large  particles, 
from  which  scattering  is  primarily  in  the  forward  direction,  the  time  variable  can 
be  eliminated  (time  would  be  equivalent  to  distance  along  the  direction  of  laser 
propagation).  For  optically  thick  media  in  which  large  numbers  of  multiple  scatters 
occur,  the  radiative  transfer  equation  reduces  to  a  diffusion  equation. 

6.4  Preliminary  Modeling 

We  have  done  preliminary  modeling  of  the  multiple  scattering  problem  by  applying 
the  radiative  transfer  equation  to  an  urban  aerosol  atmosphere,  a  medium  often 
encountered  in  our  lidar  work.  For  such  a  medium  the  extinction  coefficient  is 
typically  on  the  order  of  0.1  km-1.  Therefore,  for  lidar  paths  less  than  10  km, 
more  than  two  to  three  scatters  is  unlikely.  Furthermore,  the  particles  are  large 
compared  to  the  wavelength  of  light.  Therefore,  as  is  borne  out  by  Mie  scattering 
theory,  scattering  is  concentrated  in  a  small  cone  about  the  forward  direction. 

For  this  application  we  can  assume  that  the  transmitted  wave  propagates  in 
the  forward  direction,  attenuating  along  the  way  by  scattering  and  absorption, 
undergoes  one  backscatter,  and  propagates  back  to  the  receiver,  again  attenuating 
along  the  way  by  scattering  and  absorption.  Light  can  scatter  in  and  out  of  the 
main  propagation  path  many  times,  but  only  by  small  angles. 

For  this  application  it  is  valid  to  neglect  changes  in  polarization.  Most  scat¬ 
ters  are  small  angle,  except  for  the  one  backscatter,  which  most  likely  preserves 
polarization.  Therefore  we  may  use  Eq.  36  (only  total  specific  intensity  need  be 
considered). 

Eq.  36  can  be  reduced  from  a  6  dimensional  problem  to  a  3  dimensional  problem. 
Time  is  eliminated  as  a  degree  of  freedom;  since  the  scattering  angles  are  small 
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distance  from  the  transmitter  is  equivalent  to  time. 

There  is  azimuthal  symmetry  in  space  f  and  in  direction  s.  The  density  of  the 
medium  is  assumed  homogeneous  and  the  incident  radiation  is  assumed  azimuthally 
symmetric.  The  particles  are  assumed  to  have  an  azimuthally  symmetric  scattering 
phase  function. 

A  “solution”  of  the  radiative  transfer  equation,  for  large  particles,  is  given  by 
Ishimaru  [ll].  The  solution  is  written  as  the  following  Fourier  integral 

I(z,p,s)  =  f  j  dqe~ii4,[FQ(K,q  +  Kz)K{z,K,q)\,  (38) 

where 

F0{ictq)  =  J  dpe,K?  j  dse'‘  *I(0,p,s)  (39) 

K(z,  K,q)  =  exp  [-  ^  a(z')  (l  -  P{i j  +  k \z  -  z')))  dz']  (40) 

P(qj  =  f  dse^pW-  (41) 

In  these  expressions  p  —  xx  +  yy,  k  is  the  Fourier  transform  variable  corresponding 

to  the  variable  p,  and  q  is  the  Fourier  transform  variable  corresponding  to  the 
variable  s. 

Equation  38  indicates  that,  in  Fourier  transform  space,  the  function  K  “propa¬ 
gates”  the  initial  condition  F0  (the  Fourier  transform  of  the  initial  flux)  to  the  final 
solution.  K  is  referred  to  as  the  propagator.  Here  the  propagator  is  a  function  of 
the  Fourier  transform  P  of  the  scattering  phase  function  p. 

In  adopting  this  solution  to  our  problem  we  use  Eq.  38  to  propagate  the  incident 
radiation  at  z  =  0  to  a  height  z  =  zh.  The  result  is  multiplied  by  a  number  (the 
backscatter  coefficient).  This  result  is  then  used  as  the  initial  specific  intensity  for 
propagation  back  down  to  z  =  0.  By  multiplying  the  solution  at  z*  by  the  backscat¬ 
ter  coefficient  (instead  of  integrating  over  solid  angles)  we  assume  the  radiation  to 
specularly  reflect  with  reflectivity  given  by  the  chosen  backscatter  coefficient. 

Additional  assumptions  are 


•  The  initial  specific  intensity  is  that  of  a  Gaussian  beam: 


T/n  _  fc*Pr/0 
7(0,  p,s)  = - exp 


l  A  «})' 


(42) 
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where  p?  is  the  beam  radius  and  Or  —  v/kpr  is  rj  times  the  diffraction  limited 
transmitter  beam  spread  (l/kpr)-  is  a  free  parameter. 

•  The  scattering  phase  function  is  (Ishimaru  [11]  uses  this  phase  function) 

p(s)  =  4apWoe~a’,i  (43) 

where  W0  =  o,/at  is  the  single  particle  albedo  and  ap  gives  the  angular  spread 
of  the  single  scattered  cone  of  radiation.  ap  is  proportional  to  (D/A)2,  where 
D  is  the  diameter  of  the  particle. 


•  The  area  function  for  the  receiver  is 


P  <=  Pr 
P  >  PR 


(44) 


where  Ar  is  the  physical  area  of  the  receiver  aperture  and  pr  is  the  aperture 
radius.  The  power  collected  by  the  receiver  is  the  integral  of:  (the  specific 
intensity  at  the  receiver  multiplied  by  the  area  function). 


The  solution  to  our  problem  is  the  power  -P(z* )  collected  by  the  receiver  due  to 
scattering  at  height  z  =  zh.  Normalized  to  the  incident  power  Pt  the  solution  is 
the  following  three  dimensional  integral: 

=  (3(ZH)e-2KkaM‘  [2Wi  dA  [°°dK  [°°  dq(2ny/2  —  A)f(/c,q,  A),  (45) 

rr  Jo  Jo  Jo 

where 


f(K,q,A) 


G{q,K,A) 

Gi(9,«,A) 

Gt{q,K,A) 


Ji{q)Ji{K)  exp  --Rk1 

4 


x  exp 
x  exp 
exp 
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(q20p  -f  49rk2  +  26 p9RqK  cos  a) 


.  4 92R92 
Jo  dz\V0at{z)G{q,K, A)], 
-1 


4ap0R6p 


{Gi{q,K,  A)  +  Ga(g,  k,  A))  , 


92pq2  +  62rk 2  ^2  — — ^  +  29p6RqK  (2  -  —  j  cos  A 
62pq2  +  6rk2  ^1  -  —  ^  +  26pORqK  ^1  -  cos  A 


(46) 

(47) 

(48) 

(49) 


In  these  equations  the  following  definitions  are  used: 
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K  =  |«|, 

9=  l?1> 

A  is  the  angle  between  ic  and  q 

$T  —  v/kpT  (transmitter  beam  spread), 

Op  =  Pr/zh  (half  angle  subtended  by  receiver  at  z 
Or  —  Pt/Pr, 

Pt  is  the  transmitter  beam  radius,  and 
Ji  is  the  first  order  Bessel  function 

In  Eq.  45  the  factor  in  front  of  the  triple  integral  is  the  single  scatter  lidar  result 
multiplied  by  z\.  For  small  optical  depth  (/o*  adz  C  l)  the  triple  integral  reduces 
to  \/z\. 

Equation  45  has  been  programmed  on  an  80386  based  computer.  However,  the 
calculation  takes  several  hours  and  is  not  feasible  to  perform.  Additional  simplifi¬ 
cations  must  be  made  for  future  work. 
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A  Lidar  Experiments  at  Dripping  Spring,  NM 

A.l  Introduction 

In  June,  1988  lidar  measurements,  with  the  GL/OP  low  altitude  lidar  system,  were 
made  at  Dripping  Spring,  New  Mexico,  a  desert  region  approximately  80  miles 
south  of  Albuquerque.  The  GL/OP  low  altitude  lidar  system  measured  extinction 
coefficients  and  cloud  heights  during  this  field  program. 

Section  A. 2  discusses  the  data  analysis  techniques,  results,  and  conclusions  made 
from  these  measurements.  Section  A. 3  discusses  the  minimum  calculable  extinction 
coefficient  as  a  function  of  lidar  parameters.  We  include  this  work  here  (instead  of 
in  Section  5)  because  it  is  a  direct  result  of  the  difficulties  we  had  in  calculating  the 
small  extinction  coefficients  of  the  clear  air  of  Dripping  Spring. 

A.2  Results  and  Conclusions 

We  successfully  measured  cloud  heights  during  the  field  program.  Backscattered 
power  from  clouds  is  easily  identified  during  data  taking  either  by  the  experimenter 
or  by  the  data  acquisition  system  software. 

During  clear  air  conditions  extinction  coefficients  could  not  be  determined  be¬ 
cause  of  the  low  backscattered  power  from  the  relatively  clear  air  over  the  New 
Mexican  desert.  Instead,  we  show  that  the  measured  backscattered  power  agrees 
well  with  the  theoretically  expected  backscattered  power  from  a  purely  Rayleigh 
atmosphere.  The  total  amount  of  backscattered  power  is  due  to  Rayleigh  scattering 
off  of  atmospheric  molecules  and  to  Mie  scattering  off  of  larger  atmopsheric  aerosols. 
The  lowest  amount  of  backscatter  will  always  be  equal  to  or  greater  than  that  due 
to  Rayleigh  scattering. 

Although  the  lidar  “sees”  the  backscattered  power  from  a  Rayleigh  atmosphere, 
the  amount  is  too  low  to  be  used  to  calculate  extinction.  The  reason,  which  we 
discuss  in  the  next  section,  is  because  the  decrease  in  signal  due  to  the  “1/r*”  effect 
dominates  that  due  to  the  extinction  effect. 

A  standard  Rayleigh  atmosphere  model  was  calculated  for  each  of  the  days  of  the 
field  trip  using  measured  values  of  humidity,  pressure,  and  temperature.  The  model 
provides  values  for  the  backscatter  coefficient  0Rav  and  the  extinction  coefficient 
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dRaV  as  a  funciton  of  altitude  z.  The  backscattered  power  PRav  is  then  calculated 
in  the  following  manner: 

Pr<xV[z)  =  ^Ra^Z~  exp  (-2^  aRav(z)dz^ 

The  measured  backscattered  power  is  known  at  discrete  ranges  along  a  path  at 
some  angle  with  respect  to  the  horizontal.  Each  range  point  at  which  the  backscat¬ 
tered  power  is  measured  is  compared  to  the  calculated  backscattered  power  at  the 
same  altitude.  If  a  calculated  backscattered  power  for  a  given  measured  range 
interval  does  not  exist  then  it  is  interpolated  from  the  surrounding  points  where 
calculations  were  made.  This  is  done  up  to  an  altitude  of  15  km,  the  highest  altitude 
considered  in  the  model. 

The  lidar  system  is  not  absolutely  calibrated  and  therefore  only  relative  backscat¬ 
tered  powers  are  known.  The  theoretical  backscattered  power  is  multiplied  by  a 
constant  which  is  determined  by  linearly  fitting  it  to  the  measured  backscattered 
power  at  the  initial  range. 

The  lower  curve  labeled  “data”  in  Fig.  25a  is  a  lidar  return  with  a  cloud  at  5 
km.  The  air  below  5  km  is  clear.  The  vertical  axis  for  this  curve  corresponds  to 
a  number  proportional  to  the  measured  backscattered  power.  The  horizontal  axis 
corresponds  to  the  altitude  above  sea  level.  The  lidar  was  pointing  35°  above  the 
horizontal  for  this  return.  Dripping  Spring  is  1.8  km  above  sea  level. 

The  lower  curve  labeled  “Rayleigh  model”  in  Fig. 25b  is  the  theoretically  calcu¬ 
lated  backscattered  power  assuming  a  Rayleigh  atmosphere. 

The  upper  curve  in  Fig.  25  is  the  ratio  of  lidar  backscattered  power  to  Rayleigh 
backscattered  power.  The  vertical  axis  for  this  curve  is  just  a  dimensionless  number. 
We  see  that  the  Rayleigh  model  predicts  the  backscattered  power  fairly  well  up  to 
the  cloud.  The  cloud  produces  backscattered  power  100  times  the  level  predicted 
by  Rayleigh. 

All  other  lidar  returns  for  clear  air  conditions  show  similar  agreement  between 
data  and  Rayleigh  model. 
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Dripping  Spring,  NM:  1.8  km  above  sea  level 
lidar  elevation:  35° 


data 


Rayleigh  model 


an 


Rayleigh  model 


altitude  above  sea  level  (km) 


Figure  25:  Lower:  lidar  backscattered  power  (with  a  cloud)  and  Rayleigh  model 
backscattered  power.  Upper:  ratio  of  data  to  model. 


A. 3  Minimum  Calculable  Extinction 

A. 3.1  General  Discussion 

It  is  difficult  to  determine  small  values  of  extinction  coefficient  (a  <  0.1  km-1)  from 
lidar  data,  even  if  there  is  no  noise.  This  difficulty  arises  because  the  decrease  of  the 
lidar  return  is  nearly  the  same  as  its  decrease  due  to  the  1/r2  effect  alone  (assuming 
constant  extinction  and  backscatter  coefficient).  This  effect  is  more  severe  for  near 
ranges;  it  is  easier  to  measure  small  extinction  coefficients  at  far  ranges,  if  the  return 
signal  is  large  enough. 

The  lidar  return  P(r)  for  constant  extinction  coefficient  a  and  constant  backscat¬ 
ter  coefficient  (in  units  for  which  the  backscatter  coefficient  is  one)  is 

P(r)  =  P(r0)^exp|-2a(r-r„)|  (50) 

where  r0  is  a  reference  range  where  the  lidar  return  is  known  to  be  P(r„).  The 
decrease  of  the  lidar  return  with  the  l/r2  effect  alone  is  expressed  as 

P,*=P(r.)Tf2  (51) 

To  to  detect  the  extinction  effect  the  lidar  receiver  must  be  sensitive  to  signals 
of  at  least  the  order  of  Prt  -  P.  Figure  26  shows  a  plot  of  6  =  ( Pr»  —  P)/P(r0) 
versus  range  for  different  values  of  extinction,  with  ra  —  0.25  km.  Figure  27  shows 
6  versus  range  for  different  values  of  r0,  with  a  =  0.015  km-1. 

Figure  26  shows  that  for  Rayleigh  scattering  (a  «  0.015  km-1)  the  maximum 
value  of  6  is  approximately  2  x  10-3  times  the  power  at  the  reference  position. 
This  amount  of  power  is  difficult  to  measure  if  the  reference  power  is  in  the  range 
of  sensitivity  of  the  detector;  photomultiplier  tubes  operating  in  current  mode  are 
typically  only  sensitive  to  approximately  three  orders  of  magnitude  of  signal.  This 
amount  of  power  could  be  measured  if  the  detector  were  set  to  be  saturated  by  the 
reference  power. 

Figure  27,  however,  shows  that  the  effect  is  not  as  severe  if  the  reference  range 
is  increased.  This  is  because  the  l/r2  effect  is  determined  by  the  ratio  of  two  range 
values,  while  the  extinction  is  determined  by  the  difference  of  two  range  values.  The 
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Figure  26:  Difference  between  lidar  returns  with  and  without  decrease  due  to  ex¬ 
tinction.  Returns  are  equal  at  r0  —  0.25  km. 


range  (km) 


Figure  27:  Difference  between  lidar  returns  with  and  without  decrease  due  to  ex¬ 
tinction,  for  different  r0,  with  a  =  0.015  km-1. 


range  (km) 


Figure  28:  Difference  between  lidar  returns  with  and  without  decrease  due  to 
tinction.  Values  at  different  r0  decrease  as  1/r* 


ratio  of  two  range  values  (ri/r2,  where  r2  <  r2)  increases  as  the  range  increases,  if 
the  difference  between  the  two  ranges  remains  constant. 

Unfortunately,  moving  the  reference  range  out  may  not  help  matters  because 
the  reference  power  may  become  quite  small.  This  is  shown  in  Fig.  28,  which  is  a 
plot  of  6'  =  ( Pr2  —  P)/P„).  P0  is  the  lidar  return  at  an  initial  range  and  is  defined 
by 

P{r0)  -  P0~\ex p  [-2a  ( r0  -  r»)]  (52) 

ro 

where  r*  is  chosen  to  be  0.25  km. 

A. 3. 2  Application  to  Lidar  Data 

The  above  discussion  implies  that  we  cannot  determine  extinction  coefficients  from 
the  lidar  data  taken  during  clear  air  conditions  at  Dripping  Spring.  Previously 
we  concluded,  by  comparing  the  data  to  a  Rayleigh  model,  that  the  extinction 
coefficient  during  clear  air  conditions  is  primarily  due  to  Rayleigh  scattering  and  is 
typically  less  than  0.02  km-1. 

The  lidar  detector  was  always  set  so  that  the  return  signal  at  0.25  km  range 
was  within  the  sensitive  range  of  the  photomultiplier  tube.  The  maximum  output 
of  the  photomultiplier  tube,  in  its  linear  regime  of  operation,  is  500  mV  (into  a 
500  load).  For  clear  air  conditions  the  maximum  photomultiplier  output  (at  maxi¬ 
mum  gain)  would  be  typically  300  mV.  Figure  28  then  indicates  that  greater  than 
0.018x300  mV  =  0.54  mV  changes  in  the  lidar  return  should  have  been  detectable 
to  differentiate  between  extinction  and  1/r2  effects.  However,  a  change  of  one  bit  in 
the  analog-to-digital  converter  corresponds  to  0.49  mV  (12  bits  over  2  Volts).  Also, 
fluctuations  due  to  noise  and  background  were  typically  two  to  three  bits. 

The  problem  could  not  be  alleviated  by  increasing  the  photomultiplier  anode- 
to-cathode  voltage.  During  clear  air  conditions  this  voltage  was  always  set  to  its 
maximum  suggested  value. 

We  could  have  alleviated  the  problem  by  increasing  the  gain  of  thi.  amplifier 
which  precedes  the  analog-to-digital  converter,  thus  saturating  the  photomultiplier 
at  0.25  km.  However,  since  the  noise  fluctuations  increase  by  the  same  factor  as  the 
signal,  averaging  over  a  greater  number  of  laser  shots  would  have  been  necessary. 
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The  noise  fluctuations  decrease  linearly  with  the  square  root  of  the  number  of  laser 
shots.  Thus,  increasing  the  gain  by  a  factor  of  two  would  mean  that  four  times  as 
many  laser  pulses  would  have  to  be  averaged  over  to  get  the  same  amount  of  noise. 
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